AD-A099  337  GEOLOGICAL  SURVEY  BATON  ROUGE  LA  WATER  RESOURCES  DIV  p/6  8/8 

METHODS  AND  APPLICATIONS  OF  OI6ITAL-MODCL  SIMULATION  OF  THE  RED~CTC(U> 
MAY  80  A  H  LUDWIG*  J  E  TERRY 

UNCLASSIFIED  USGS/WRD/WRI-8X/037  NL 


„ATION.  RED  RIVER  ALLUVIAL  A 


CO 

CO 

o> 


METHODS  AND  APPLICATIONS  OF  DIGITAL-MODEL  SIMULATION  OF  THE 


RED  RIVER  ALLUVIAL  AQUIFER,  SHREVEPORT  TO  THE  MOUTH  OF  THE 
BLACK  RIVER,  LOUISIANA 


U.S.  GEOLOGICAL  SURVEY 


Water-Resources  Investigations  79-114 


Prepared  in  cooperation  with  the 


U.S.  Army  Corps  of  Engineers 
and  the 

U.S.  Soil  Conservation  Service 


£ 

M 


81  5 


26  024 


30371-10) 


REPORT  DOCUMENTATION 
PAGE 


J.  REPORT  NO. 


IJSGS/1\'R])/1\KI  -  81/037 


/ 


4.  Tttl.  and  Subtitle 

METHODS  AND  APPLICATIONS  OF  DIGITAL-MODEL  SIMULATION  OF  THE 
RED  RIVER  ALLUVIAL  AQUIFER,  SHREVEPORT  TO  THE  MOUTH  OF  THE 
BLACK  RIVER,  LOUISIANA  _ 


7.  Authors) 

A.  H.  Ludwig  and  J.  E.  Terry 

9.  Performing  Organization  Name  and  Address 


U.S.  Geological  Survey,  Water  Resources  Division 
P.O.  Box  66492 

Baton  Rouge,  Louisiana  70896 


12.  Sponsoring  Organization  Name  and  Address 

U.S.  Army  Corps  of  Engineers 
New  Orleans  District 
P.O.  Box  60267 

New  Orleans,  Louisiana  70160 


1.  Recipient's  Acc.t.ion  No. 


5.  Report  Dot* 

May  1980 

6. 

8.  Performing  Organization  Rapt.  No. 

USGS/WRI-79- 114  V 

10.  Project  /Task  /  Work  Unit  No. 

11.  Contract(C)  or  Grant(G)  No. 

(C) 

(G) 


13.  Type  of  Report  &  Period  Covered 

Final 


14. 


15.  Supplementary  Notes 


Prepared  in  cooperation  with  the  U.S.  Army  Corps  of  Engineers  and  the  U.S.  Soil 
Conservation  Service 


16.  Abstract  (Limit:  200  words.) 

. .  The  Red  River  Waterways  Project  provides  for  the  construction  of  five  locks  and 

dams  on  the  Red  River  from  the  Mississippi  River  to  Shreveport,  La.  The  methodology 
used  by  the  U.S.  Geological  Survey  in  studying  the  effects  of  the  navigation  pools  on 
the  ground-water-f low  regime  involved  digital  modeling  of  steady-  and  nonsteady-state 
conditions.  The  steady-state  model,  GWFLOW,  computes  the  head  response  in  an  aquifer 
due  to  various  boundary  conditions.  The  nonsteady-state  model,  SUPERMOCK,  was  de¬ 
signed  to  simulate  transient  stress  and  response  in  an  alluvial-flow  system. 

In  addition  to  the  simulation  models,  several  computer  programs  were  developed 
during  the  study  to  aid  in  the  preparation  of  field  data  for  input  to  the  models  and 
in  the  calibration  of  the  models.  Calibration  techniques  unique  to  each  of  the 
models  were  developed  for  the  investigation. 


17.  Document  Analysis  a.  Descriptors 

*Groundwater ,  *Aquifer  characteristics,  ^Alluvial  aquifers.  Groundwater  movement. 
Model  studies. 


b.  Identifiers /Open-Ended  Terms 

*Red  River  Valley,  Digital  modeling. 

Approved  for  public  release;  distribution 
unlimited. 

18.  Availability  Statement 


Release  unlimited 


19.  Security  Class  (This  Report) 

UNCLASSIFIED 

30.  Security  Cl*u  (Thi»  P»*»> 

UNCLASSIFIED 


31.  Nu  of  P»«*« 

no 

33.  Price 


(So*  ANSt-Z3S.lt) 


Sn  Instructions  on  Ro.orio 


'/ 


I  373  («-»n 
(formerly  NTIS-35) 
Department  of  Commerce 


METHODS  AND  ^APPLICATIONS  OF  DIGITAL-MODEL_SIMULATION  OF  THE 
RED  RIVER  ALLUVIAL  AQUIFER,  SHREVEPORT  TO  THE  MOUTH  OF  THE 
BLACK  RIVER,  LOUISIANA  . 


By  A.  H.  Ludwig  and  J.  E.  Terry 


i)  Pi 


US.  GEOLOGICAL  SURVEY 


^  /ty  qj  j  Q0 


Water-Resources  Investigations  79-114 


Prepared  in  cooperation  with  the 

U  S  Army  Corps  o>  Engineers 


and  the 

U  S  Soil  Conservation  Service 


/f  /WKlj 

T 


I- 

n  i-'xi* 


May  1980 


UNITED  STATES  DEPARTMENT  OF  THE  INTERIOR 


CECIL  D.  ANDRUS,  Secretary 
GEOLOGICAL  SURVEY 
H.  William  Menard,  Director 


For  additional  information  write  lo 

U.C  (jeological  Survey 
P  O.  Box  6649? 

Baton  Rouge.  Louisiana  70896 


J 


CONTENTS 


Factors  for  converting  inch-pound  units  to  International  System  (SI) 

of  metric  units - 

Abstract - 

Introduction - 

Background  of  the  investigation - 

Purpose  and  scope - 

Data  collection - 

Modeling  the  hydrologic  system - 

Conceptual  model - 

Digital  model - 

Nonsteady  state - 

Steady  -state - 

Preparation  of  digital-model  input - 

Nonsteady-state  model - 

Root  depth - 

Land-surface  elevation - 

Initial  potentiometric  surface  and  water  table - 

Observed  potentiometric  surface  and  water-table  elevations 

Transmissivity - 

Conductivity  of  the  upper  confining  layer - 

Relation  of  evapotranspiration  to  depth  to  water - 

Thickness  of  streambed  material - 

Specific  yield  and  storage  coefficient - 

Precipitation  and  potential  evapotranspiration - 

River  stage - 

Steady-state  model - 

Transmissivity - 

Change  in  evapotranspiration  with  change  in  potentiometric 

surface - 

Thickness  of  streambed  material - 

Head  conditions  in  confined  aquifer  (node-level  map) - 

Changes  in  stream  stage - 

Calibration  and  verification  of  the  nonsteady-state  model - 

Model  output - 

Continuing  studies - 

Selected  references - 

At  tachments - 

Attachment  A,  AVERAGE  program - 

Attachment  B,  ATMOFLUX  program - 

Attachment  C,  POTEET  program - 

Attachment  D,  RIVCHANGE  program - 

Attachment  E,  TRIBCHANGE  program - 

Attachment  F,  DELETDELH  program - 


ILLUSTRATIONS 


Page 


igure  1.  Map  showing  planned  navigation  features.  Red  River  Waterways 

Project -  3 

2.  Map  showing  location  of  project  area  and  model  coverage -  4 

3A-E.  Maps  showing  the  data-collection  network: 

A.  Lock  and  Dam  1  area -  7 

B.  Lock  and  Dam  2  area -  8 

C.  Lock  and  Dam  3  area -  9 

D.  Lock  and  Dam  4  area -  10 

E.  Lock  and  Dam  5  area -  11 

4.  Idealized  hydrogeologic  section  of  the  Red  River  Valley -  14 

5.  Diagram  showing  relation  between  soil-moisture  accounting, 

vertical-flow,  and  horizontal-flow  components  of  SUPERMOCK 
program —  -  17 

6.  Flow  diagram  of  digital-model  procedure  for  nonsteady-state 

analysis - 23 

7.  Example  of  alphameric  root-depth  map -  25 

8.  Map  showing  topographic  coverage  of  the  project  area -  26 

9.  Example  of  alphameric  transmissivity  map -  29 

10.  Trilinear  graph  of  soil-classification  scheme  showing 

hydraulic-conductivity  values  for  soil  classes  used  in 
ATMOFLUX  program -  31 

11.  Example  of  alphameric  map  of  vertical  hydraulic  conductivity 

and  explanation  of  symbols -  32 

12.  Example  of  GWETO  functions  for  computation  of  evapotran- 

spiration -  34 

13.  Example  of  alphameric  streambed-thickness  map  for  nonsteady- 

state  analysis -  36 

14.  Example  of  alphameric  specific-yield  and  storage-coefficient 

map  and  explanation  of  symbols -  38 

15.  Flow  diagram  of  digital-model  procedure  for  steady-state 

analysis -  41 

16.  Example  of  alphameric  map  of  AET/AH  and  explanation  of 

symbols -  44 

17.  Example  of  alphameric  streambed-thickness  map  for  steady- 

state  analysis -  45 

18.  Example  of  node-level  map -  46 

19.  Plan-and-prof ile  views  of  a  segment  of  river  channel  showing 

the  method  for  computing  stage  change -  48 

20.  Examples  of  spring  and  fall  calibration  charts  from  DATE 

program -  50 

21.  Example  of  accretion-summation  chart -  53 

22.  Example  of  computed  output  from  nonsteady-state  model -  54 

23.  Map  showing  average  preconstruction  potent iometric  surface. 

Lock  and  Dam  3  area -  56 

24.  Example  of  computed  output  from  steady-state  model  showing 

a  section  of  head-change  map -  57 

25.  Contour  map  showing  computed  head  change.  Lock  and  Dam  3 

area -  58 


IV 


Page 


Figure  26.  Map  showing  computed  average  postconstruction  potentiometric 

surface,  Lock  and  Dam  3  area -  59 

27.  Generalized  chart  showing  relationship  of  digital  programs 

that  prepared  data  for  input  to  SUPERMOCK  and  GWFLOW 

models -  62 

28.  Example  of  output  from  AVERAGE  program -  72 

29.  Example  of  output .Jrrom  DELETDELH  program -  103 


TABLES 


Page 

Table  1.  Specifications  for  lock  and  dam  arrangements  studied  in  the 

investigation -  16 

2.  Input  data  for  AVERAGE  program -  64 

3.  AVERAGE  program  listing -  66 

4.  Input  data  for  ATMOFLUX  program -  74 

5.  ATMOFLUX  program  listing -  76 

6.  Input  data  for  POTEET  program -  80 

7.  POTEET  program  listing - 82 

8.  Input  data  for  RIVCHANGE  program -  86 

9.  RIVCHANGE  program  listing -  88 

10.  Input  data  for  TRIBCHANGE  program -  94 

11.  TRIBCHANGE  program  listing -  96 

12.  Input  data  for  DELETDELH  program -  100 

13.  DELETDELH  program  listing - 101 


V 


FACTORS  FOR  CONVERTING  INCH- POUND  UNITS  TO  INTERNATIONAL  SYSTEM  (SI)  OF 

METRIC  UNITS 

A  dual  system  of  measurements — inch-pound  units  and  the  International 
System  (SI)  of  metric  units — is  given  in  this  report.  SI  is  a  consistent 
system  of  units  adopted  by  the  Eleventh  General  Conference  of  Weights  and 
Measures  in  1960.  The  conversion  factors  for  terms  used  in  this  report  are 
as  follows: 


Multiply  inch-pound  unit  By 

acre  4,047 

inch  (in.)  25.40 

inch  per  day  (in/d)  25.40 

foot  (ft)  0.3048 

foot  per  day  (ft/d)  0.3048 

foot  per  year  (ft/yr)  0.3048 

foot  squared  per  day  (ft2/d)  0.09290 

mile  (mi)  1.609 

square  mile  (mi  )  2.590 


To  obtain  SI  unit 
square  meter  (m^) 
millimeter  (mm) 
millimeter  per  day  (mm/d) 
meter  (m) 

meter  per  day  (m/d) 
meter  per  year  (m/vear) 
meter  squared  per  day  (m-/d) 
kilometer  (km) 
square  kilometer  (km  ) 
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METHODS  AND  APPLICATIONS  OF  DIGITAL-MODEL  SIMULATION  OF  THE 
RED  RIVER  ALLUVIAL  AQUIFER,  SHREVEPORT  TO  THE 
MOUTH  OF  THE  BLACK  RIVER,  LOUISIANA 


By  A.  H.  Ludwig  and  J.  E.  Terry 


ABSTRACT 

The  Red  River  Waterways  Project  of  the  U.S.  Army  Corps  of  Engineers 
provides  for  the  construction  of  a  series  of  locks  and  dams  on  the  Red  River 
from  the  Mississippi  River  to  Shreveport,  La.  The  locks  and  dams  will  cause  a 
permanent  rise  in  the  level  of  the  river,  creating  changes  in  the  ground-water 
flow  system.  The  U.S.  Geological  Survey,  in  cooperation  with  the  Corps  and 
the  U.S.  Soil  Conservation  Service,  began  an  investigation  in  1968  to  studv 
the  effects  of  the  planned  navigation  pools  on  the  ground-water  flow  regime. 

The  Red  River  downstream  from  Shreveport  flows  through  an  alluvial  valley 
that  ranges  from  2  to  12  miles  (3.2  to  19  kilometers)  in  width.  Along  the 
thalweg  of  the  valley,  the  alluvium  ranges  from  75  to  200  feet  (23  to  61 
meters)  in  thickness  and  is  composed  of  a  silt  and  clay  layer,  underlain  by  a 
coarse  sand  and  gravel  aquifer.  The  aquifer  is  hydraulically  connected  in 
varying  degrees  to  the  Red  River  and  its  major  tributaries. 

The  methods  used  in  the  investigation  involved  digital  modeling  of  steadv- 
and  nonsteady-state  conditions.  The  nonsteady-state  model,  utilizing  a  program 
called  SUPERMOCK,  was  designed  to  simulate  transient  stress  and  response  in  a 
ground-water  flow  system  that  includes  a  water  table  in  a  confining  laver 
above  an  artesian  aquifer.  The  steady-state  model,  utilizing  a  program  called 
GWFLOW,  computes  the  head  response  in  an  aquifer  due  to  various  boundary 
conditions. 

Principal  data  requirements  for  the  models  include  climatic  data,  defini¬ 
tion  of  the  hydraulic  characteristics  of  tiie  upper  confining  layer  and  aquifer, 
water-table  levels  in  the  upper  confining  layer  and  potontiometric  levels  in 
the  aquifer,  and  stream-stage  data  for  the  Red  River  and  its  tributaries. 

In  addition  to  the  simulation  models,  several  computer  programs  were 
developed  to  aid  in  preparation  of  data  and  in  the  calibration  of  the  models. 
The  programs  were  designed  to  compute  tiie  harmonic-mean  water  level  at  each 
observation  well  (AVERAGE),  compute  the  harmonic-mean  conductivity  for  layered 
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materials  and  the  potential  upward  movement  of  water  due  to  evapotranspiration 
at  the  land  surface  (ATMOFLUX),  compute  daily  evapotranspiration  (POTEET), 
provide  main-stem  and  tributary  stream-stage  data  sets  for  the  nonsteady-state 
model  (RIVCHANCE  and  TRIBCHANCE),  and  to  compute  the  change  in  the  rate  of 
evapotranspiration  due  to  a  change  in  protentiometric  head  (DELETDELH). 

Calibration  techniques  unique  to  each  of  the  models  were  developed  for 
the  investigation.  The  calibration  procedure  for  the  nonsteady-state  model 
involved  reproducing,  by  manipulation  of  model  parameters  within  plausible 
limits,  observed  water-table  and  potentiometr ic  levels  while  maintaining 
reasonable  limits  on  the  rate  of  accretion  to  the  aquifer. 


INTRODUCTION 


Background  of  the  Investigation 


The  Red  River  Waterways  Project  of  the  U.S.  Army  Corps  of  Engineers  was 
authorized  by  the  90th  Congress  in  the  Rivers  and  Harbors  Act  of  1958.  Project 
plans  include  a  9-  by  200-foot  (2.7-  by  61-m)  navigation  channel,  beginning  at 
the  confluence  of  the  Red  and  Mississippi  Rivers  and  winding  northwestward 
along  the  present  course  of  the  Red  River  to  Shreveport,  La.  From  Shreveport 
the  channel  will  follow  Twelvemile  and  Cypress  Bayous  to  a  point  in  Lake  O' 
the  Pines  Reservoir  near  Daingerf ield ,  Tex.  (fig.  1).  A  series  of  eight  locks 
and  dams  will  be  required  to  provide  the  navigation  depths  and  the  necessary 
225-foot  (69-m)  lift  from  the  Mississippi  River  to  the  head  of  navigation. 


The  natural  ground-water  flow  system  in  the  Red  River  alluvial  valley 
will  be  altered  by  the  formation  of  navigation  pools  except  at  locks  7  and  9, 
which  are  to  be  built  into  existing  dams  on  Caddo  Lake  and  Lake  O'  the  Pines. 
Predominant  effects  of  the  navigation  pools  on  the  ground-water  regime  will  be 
a  rise  in  water  levels  and  changes  in  the  ground-water  flow  pattern.  In  April 
1963,  at  the  request  of  the  Corps  of  Engineers,  the  U.S.  Geological  Survey 
began  a  preliminary  study  of  the  preconstruction  and  postconstruction  ground- 
water  conditions.  The  study  characterized,  using  available  data,  the  existing 
ground-water  conditions  in  the  valley  and  provided  steady-state  projections  of 
the  effects  of  proposed  navigation  structures  on  ground-water  levels.  The 
projections  were  made  with  the  aid  of  an  analog  model. 


In  1968  tlie  Corps  requested  that  the  Geological  Survey  refine  the  projec¬ 
tions  made  in  the  earlier  study  and  that  a  continuing  ground-water  data- 
collection  program  in  the  Red  River  Valley  be  established.  The  study  area  was 
the  alluvial  valley  from  the  confluence  of  the  Red  and  Black  Rivers  to  Shreve¬ 
port,  La.,  a  distance  of  241  river  miles  or  388  km  (fig.  2).  The  Corps  of 
Engineers  considered  several  arrangements  of  either  five  or  six  locks  and  dams 
within  this  reach  of  the  river.  An  arrangement  of  five  locks  and  dams,  known 
as  tlie  B- 3  modified  plan,  was  considered  the  most  feasible  plan  of  construction 


The  effects  of  increased  river  stages,  caused  by  the  formation  of  naviga¬ 
tion  pools,  on  the  ground-water  regime  were  projected  for  steady-  and  nonsteadv 
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Planned  navigation  features,  Red  River  Waterways  Project. 


state  conditions.  The  steady-state  change  in  potentiometric  surface  was 
projected  from  the  average  postconstruction  river  stages.  Nonsteady-state 
projections,  used  to  determine  the  effect  of  the  project  on  agriculture, 
were  made  for  specific  periods  within  a  typical  calendar  year.  Because  of 
the  size  of  the  project  area  and  the  complexity  of  the  flow  regime,  digital¬ 
modeling  techniques  were  used. 

Results  of  the  steady-  and  nonsteady-state  analyses  for  each  of  the 
five  lock-and-dam  areas  were  provided  to  the  Corps  of  Engineers,  1975-76,  in 
a  series  of  five  administrative  reports  that  were  later  released  to  the  open 
file  (Ludwig,  1979a,  b;  Ludwig  and  Reed,  1979;  Ludwig  and  Terry,  1979a,  b). 
Average  depths  to  the  water  table  for  specific  periods  of  interest  were 
prepared  for  the  U.S.  Soil  Conservation  Service  on  punched  IBM  computer 
cards.  Basic-data  reports  containing  ground-water  quality  analyses  (Ludwig, 
1974)  and  ground-water  levels  through  June  1975  (Stephens,  1976)  were  pub¬ 
lished  as  open-file  reports. 


Purpose  and  Scope 

The  purpose  of  this  report  is  to  describe  the  methods  used  in  the  study 
and  to  show  their  application  to  the  Red  River  Waterways  Project.  The  dis¬ 
cussion  is  intended  to  be  sufficiently  detailed  that  the  reader  can  obtain  a 
basic  understanding  of  the  methodology  employed  in  the  study.  The  discussion 
covers  (1)  development  and  management  of  the  basic-data  network  and  the  types 
of  data  collected,  (2)  conceptualization  of  the  geohydrology  of  the  study 
area,  (3)  descriptions  of  predictive  models  used  and  data  requirements  of 
the  models,  (4)  presentation  of  peripheral  digi tal-computer  programs  used  to 
generate  or  manipulate  data  for  use  in  the  models,  (5)  calibration  of  the 
models,  (6)  descriptions  of  output  from  the  models,  and  (7)  possible  utiliza¬ 
tion  of  the  calibrated  Red  River  models  for  other  uses.  Examples  of  program 
input  and  output  (taken  from  analyses  of  Lock  and  Dam  3  area)  are  shown. 


DATA  COLLECTION 

The  objective  of  the  data-collection  program  for  the  Red  River  study  was 
to  obtain  the  data  necessary  for  the  determination  of  the  hydrologic  charac- 
terics  of  the  flow  regime  in  the  Red  River  alluvium  and  the  climatic  factors 
and  agricultural  practices  which  affect  it.  To  accomplish  these  objectives, 
work  activities  were  divided  among  the  participating  agencies  as  follows:  The 
Geological  Survey  mapped  the  principal  hydrologic  boundaries,  inventoried 
existing  wells  suitable  for  periodic  measurements,  drilled  test  holes  and 
installed  observation  wells,  analyzed  samples  of  alluvial  material  for 
hydraulic  conductivity  and  grain  size,  installed  and  operated  a  series  of 
surface-water  gages  on  tributary  streams,  and  analyzed  ground-water  samples 
from  selected  wells  for  chemical  constituents.  The  U.S.  Soil  Conservation 
Service  installed  shallow  piezometers  at  observation-well  sites,  monitored 
crop-observation  plots  to  establish  the  relationship  between  yield  and  soil- 
moisture  conditions,  mapped  soil  profiles,  inventoried  land-use  practices, 
and  measured  water  levels  in  the  network  of  Geological  Survey  and  Soil 
Conservation  Service  observation  wells  and  piezometers.  The  Corps  of  Engi- 


neers  provided  average  preconstruction  and  postconstruction  stage  profiles 
of  the  Red  River  to  be  used  in  developing  input  to  the  steady-state  model. 

The  Corps  also  provided  time-variant  preconstruction  and  postconstruction 
stage  data  in  the  form  of  5-day  averages  at  2-mile  (3.2-km)  increments  for 
the  period  December  1967  to  September  1973  for  the  entire  reach  of  the  Red 
River  in  the  project  area. 

The  test-drilling  program  conducted  by  the  Geological  Survey  was  com¬ 
pleted  during  a  series  of  field  sessions  from  1968  to  1971.  Approximately  350 
test  holes  were  drilled  in  the  valley,  from  Shreveport  to  the  mouth  of  the 
Black  River.  Test  holes  were  drilled  with  solid-stem  power-auger  drilling 
equipment,  and  soil  samples  were  collected  at  selected  depths  for  analyses  of 
hydraulic  conductivity  and  particle-size  distribution.  Most  of  the  test  holes 
were  drilled  and  logged  through  the  entire  alluvial  section  and  into  the 
underlying  Tertiary  bedrock.  The  test  holes  were  cased  with  l'i-inch  (32-mm) 
galvanized-iron  pipe  and  screened  with  3-foot  (0.9-m),  60-gage  well  screens. 

The  screens  were  set  opposite  coarse  sand  and  gravel  at  depths  ranging  from  20 
to  140  ft  (6  to  43  m)  below  the  land  surface.  The  locations  of  the  observation 
wells  are  shown  in  figures  3A-E. 

In  the  vicinity  of  the  proposed  construction  sites  and  along  the  river, 
the  wells  are  more  closely  spaced  in  anticipation  of  greater  variations  in 
water  levels  in  these  areas.  At  greater  distances  from  the  river,  fewer  wells 
are  required.  The  amount  of  pumpage  from  the  alluvium  is  small;  therefore, 
where  little  change  was  expected,  the  data  from  a  particular  well  could  be 
extrapolated  over  a  relatively  large  area.  The  density  of  wells  ranged  from 
one  well  per  square  mile  (2.6  km2)  in  the  vicinity  of  the  locks  and  dams  to 
about  one  well  per  3  mi2  (7.8  km2)  elsewhere  in  the  valley. 

Shallow  piezometers  were  placed  adjacent  to  most  of  the  observation  wells 
to  obtain  data  on  the  position  of  the  water  table  in  the  upper  confining 
layer.  The  piezometers  consisted  of  lengths  of  3/4-inch  (19-mm)  galvanized- 
iron  pipe,  driven  into  the  ground  to  selected  depths  ranging  from  1  to  20  ft 
(0.3  to  6.1  m)  below  the  land  surface.  The  lower  end  of  the  pipe  was  left 
open  to  the  soil  to  allow  movement  of  water  into  and  out  of  the  pipe.  Two  to 
five  piezometers  were  installed  at  each  observation-well  location,  depending 
on  the  variations  in  lithology  in  the  upper  section. 

Water-level  measurements  in  all  observation  wells  and  piezometer  tubes 
were  made  monthly  by  Soil  Conservation  Service  personnel.  Digital  recorders 
were  installed  on  16  wells  in  the  study  area.  Fourteen  of  the  wells  were  near 
the  Red  River  to  provide  daily  water-level  data  for  the  computation  of  aquifer 
diffusivity.  In  addition,  water  samples  were  collected  from  all  of  the 
observation  wells  at  the  time  of  installation  and  from  many  piezometer  tubes 
and  analyzed  for  chemical  quality. 

Stream-stage  data  were  collected  from  a  network  of  45  continuous  recorders, 
staff  gages,  and  wire-weight  gages  (figs.  3A-E).  Most  of  the  gages  were  part 
of  the  regular  surface-water  data-collection  network  operated  by  the  Geological 
Survey  and  the  Corps  of  Engineers.  However,  14  additional  gages  were  installed 
at  intervals  along  tributary  streams  between  existing  recording  gages  and  on 


9m9  station  (continuous  -  record)  A  Stream  -  stage  site  (nonrecordmj) 


Stream- gaging  station  (continuous 


lakes  in  the  valley.  These  gages  provide  supplementary  data  for  the  determina¬ 
tion  of  stream  profiles.  Climatic  data,  including  maximum  and  minimum  dailv 
temperature  and  daily  precipitation,  were  obtained  from  five  National  Weather 
Service  stations  in  and  near  the  valley  (figs.  3A-E). 


MODELING  THE  HYDROLOGIC  SYSTEM 


Conceptual  Model 


The  Red  River  flows  southeastward  through  central  and  northwestern  Loui¬ 
siana.  From  Shreveport  to  the  vicinity  of  Marksville,  the  river  is  confined  in 
a  valley  ranging  from  2  to  12  mi  (3.2  to  19  km)  in  width.  The  uplands  bordering 
tlu-  valley  rise  as  much  as  150  ft  (46  m)  above  the  general  level  of  the  valley. 
Downstream  from  Marksville,  the  Red  River  Valley  and  the  Mississippi  River 
vallev  merge  to  form  the  broad  Mississippi  River  alluvial  plain.  The  flood 
plain  is  characterized  by  very  low  relief,  meandering  stream  courses,  oxbow 
lakes,  and  other  alluvial  features.  The  dominant  features  are  natural  levees, 
which  form  the  topographic  highs,  and  backswamps,  which  are  the  topographic 
lows.  Tlie  natural  levees  rise  from  10  to  20  ft  (3  to  6  m)  above  the  adjoining 
backswamps.  Natural  levees  occur  along  abandoned  channels  of  the  Red  River 
and  on  tributary  streams,  as  well  as  along  the  present  course  of  the  river. 


Elevations  in  the  valley  range  from  40  f t  (12.2  m)  above  moan  sea  level 
(now  generally  referred  to  as  National  Geodetic  Vertical  Datum  of  1929),  near 
the  confluence  of  the  Red  and  Black  Rivers,  to  170  ft  (52  m)  above  sea  level, 
at  Shreveport. 


The  average  annual  precipitation  in  the  valley  ranges  from  57  in.  (1,448 
mm),  at  Alexandria,  to  43  in.  (1,092  mm),  at  Shreveport.  The  greatest  precipi¬ 
tation  generally  occurs  in  April  and  May,  and  the  least  in  September  and 
October.  The  climate  of  the  area  is  classified  as  humid;  that  is,  precipitation 
equals  or  exceeds  potential  evapotransp  i  rat  ion.  Favorable  climatic  conditions 
and  rich  soil  support  abundant  vegetal  growth.  In  general,  row  crops,  prin¬ 
cipally  cotton  and  soybeans,  are  grown  on  the  natural  levees.  The  lower 
levels  of  the  natural  levees  arc  used  mainlv  for  pasture  or  soybeans,  and  t  lie 
backswamp  areas  are  mostly  forested. 


Formations  of  Tertiarv  age  underlie  the  vallev  alluvium  and  crop  out 
along  the  vallev  walls.  The  beds  are  composed  primarily  oi  clav,  hut  locally 
thev  contain  sand  lenses.  The  beds  form  a  nearly  impermeable  boundary  to  the 
alluvial  aquifer.  In  many  places.  Pleistocene  terrace  deposits  overlie  the 
Tertiary  deposits  in  the  upland.  The  terrace  deposits,  which  arc  remnants  o! 
older  and  higher  flood  plains  ot  the  Red  River,  are  most  prevalent  in  the 
lowe r  end  ol  the  vallev,  where  thev  are  as  much  .is  200  tt  (*>1  m)  thick.  The 
Marksville  Prairie  is  a  terrace  remnant  in  the  Red  River  flood  plain.  Phe 
terrace  deposits  are  composed  of  a  heterogeneous  sequence  ot  sand,  silt,  and 
clay.  Gravel  lavers  occur  in  the  terrace  deposits  and  local Iv  are  the  soon  e 
of  large  quantities  of  water. 
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The  alluvium  in  the  valley  generally  ranges  from  about  75  ft  (23  m)  in 
thickness,  in  the  upper  end  of  the  area,  to  about  200  ft  (61  m),  downstream 
from  Mark.sville.  The  alluvium  can  be  divided  into  two  segments:  a  lower  unit 
or  aquifer,  which  is  composed  of  coarse  sand  and  gravel  grading  upward  to  fine 
sand,  and  an  upper  confining  layer,  which  is  composed  of  clay,  silt,  and  fine 
sand.  The  upper  confining  layer  averages  about  30  ft  (9.1  m)  in  thickness  and 
ranges  from  a  few  feet  to  140  ft  (43  m).  The  aquifer  ranges  from  5  ft  (1.5  m) 
in  thickness  beneath  some  channel-fill  and  backswamp  deposits  to  150  ft  (46  m) 
in  the  lower  end  of  the  valley.  The  thicknesses  of  the  two  segments  vary  from 
place  to  place.  Differences  of  as  much  as  100  ft  (30  m)  in  the  thickness  of 
the  upper  confining  layer  within  short  distances  have  been  noted  in  Lock  and 
Dam  1  area.  To  a  lesser  extent,  variations  in  thickness  occur  at  many  places 
in  the  valley,  primarily  as  the  result  of  fine-grained  deposition  in  former 
channels  of  the  Red  River. 

Throughout  the  Red  River  Valley,  the  Red  River  and  its  major  tributaries 
are  hvdraul icallv  connected  in  varying  degrees  to  the  Red  River  alluvial 
aquifer.  Therefore,  changes  in  stream  stages  resulting  from  the  construction 
of  the  proposed  locks  and  dams  would  induce  similar  changes  in  the  potent io- 
metric  surface  of  the  aquifer.  The  potentiometric  surface  refers  to  the  level 
to  which  water  will  rise  in  wells  tapping  the  aquifer.  Also,  throughout  the 
Red  River  Valley  a  water  table  exists  as  the  upper  surface  of  the  zone  of 
saturation  in  the  fine-grained  material  above  the  aquifer.  The  altitude  of 
the  water  table  at  any  point  is  a  function  of  the  transient  flow  through  the 
fine-grained  material  above  the  aquifer  and  the  transient  head  in  the  aquifer. 
Therefore,  induced  changes  in  the  position  of  the  potentiometric  surface  would 
indirectly  cause  changes  in  the  position  of  the  water  table. 

Rainfall  on  the  flood  plain  is  the  primary  source  of  recharge  for  the 
alluvial  aquifer.  Moisture  reaches  the  aquifer  indirectly  bv  infiltrating  the 
fine-grained  material  in  the  confining  layer  above  the  aquifer.  An  unknown, 
but  probably  very  small,  amount  of  recharge  is  derived  from  the  formations  of 
Tertiarv  age  that  underlie  and  flank  the  valley.  Most  of  the  water  moving 
downgradient  through  the  terrace  deposits  is  discharged  into  the  tributary 
streams  that  flow  along  the  margin  of  the  valley. 

Water  levels  in  most  wells  tapping  the  aquifer  rise  above  the  base  of 
tile  fine-grained  material  overlving  the  aquifer,  an  indication  that  the  water 
is  under  confined  or  semiconfined  conditions.  A  zone  of  saturation  in  the 
upper  fine-grained  material,  extending  from  near  the  land  surface  down  to  the 
aquifer,  indicates  the  presence  of  water-table  conditions.  These  two  con¬ 
ditions  i'xist  simultaneously  because  of  the  great  difference  in  hvdraul ic 
conductivity  between  tile  fine-grained  material  overlving  the  aquifer  and  the 
aquifer  itself.  The  position  of  the  water  table  may  he  either  above  or  below 
the  potentiometric  level  in  the  aquifer,  as  reflected  bv  the  direction  of  the 
resultant  vertical  flow  in  the  fine-grained  material  between  the  water  table 
and  the  top  of  t  tie  aquifer.  Accretion,  as  del  ined  bv  Stallman  ( 1 V  5h ) ,  is  the 
rati'  at  which  water  is  gained  or  lost  through  the  aquifer  surface  in  response 
to  precipitation  and  evapot ranspi rat  ion.  Positive  accretion  or  recharge  takes 
place  where  the  vertical  hvdraul ic  gradient  is  downward.  Conversely,  negative 
accretion  or  discharge  takes  place  where  the  vertical  hvdraul ic  gradient  is 
upward . 
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The  natural  movement  of  water  in  the  alluvium  is  toward  discharge  points 
along  the  Red  River  and  its  tributaries  in  the  valley.  Because  pumpage  of 
water  from  wells  is  not  significant,  water  levels  in  the  alluvium  fluctuate  in 
response  to  seasonal  variations  in  precipitation,  evapo t ranspi rat i on ,  and  to 
changes  in  river  stage. 

The  recharge,  movement,  and  discharge  of  water  from  the  alluvial  aquifer 
are  shown  graphically  in  the  idealized  section  in  figure  4.  The  direction  of 
water  movement,  indicated  by  arrows,  stiows  that  the  aquifer  is  being  recharged 
in  zone  1  where  the  gradient  is  downward  through  the  clay  and  silt.  Discharge 
takes  place  to  the  Red  River  and  vertically  upward  in  zone  2.  The  flow  con¬ 
ditions  shown  in  the  diagram  may  change.  At  any  given  location,  the  rate  of 
accretion  is  neither  constant  nor  in  the  same  direction  at  all  times.  Seasonal 
weather  changes,  changes  in  river  stage,  and  pumping  may  cause  variations  in 
the  magnitude  and  direction  of  water  movement  in  the  aquifer. 
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Figure  4  Idouli/od  hydrogeologic  section  of  the  Red  Riw'r  V.illoy 


Two  typos  of  digital  models  woro  used  in  tin1  nnalvsos.  A  s Load v-st a t o 
modol,  liWKLOW  (Bedinger  and  othors,  1973),  was  usod  to  provide  projections  of 
changes  In  the  potent iomet r i c  surface.  A  nonsteady-state  model,  SIM’KKMOCK 
(Reed  and  others,  197b),  was  used  to  simulate  fluctuations  of  the  head  in  the 
aquifer  and  the  water  table.  For  purposes  of  analysis,  the  project  area  was 
divided  into  I ive  overlapping  model  areas.  Each  area  contained  one  or  more  of 
the  proposed  loek-and-dam  construction  sites.  These  areas  are  identified  bv 
referring  to  a  particular  lock-and-dam  area  (I ig.  2).  To  aid  the  Corps  of 
Engineers  in  determining  the  best  arrangement  of  locks  and  dams,  steadv-state 
analyses  were  run  for  all  alternate  plans,  including  the  B- 3  modified  plan. 
Specifications  for  dam  locations  and  pool  elevations  for  the  plans  considered 
are  shown  in  table  1.  The  nonsteadv-state  model  was  used  to  make  projections 
for  the  15-1  modified  plan  onlv. 

The  framework  lor  the  digital  models  consisted  of  a  rectangular  grid  of 
34  rows  and  >30  columns  superimposed  on  a  map  of  the  area  having  a  scale  of 
1:62, 300.  Tile  spacing  between  each  intersection  (node)  in  the  .grid  represented 
a  distance  of  0.3  mi  (0.8  km).  Thus,  each  modol  represented  a  17-  bv  40-mile 
(27-  bv  (>4-km)  area.  Five  such  models  were  vised,  each  representing  a  lock- 
and-dam  area,  to  cover  the  190-mile  (30b- km)  reach  of  navigation  channel  in 
the  studv  area  (figs.  2,  3A-E). 

The  examples  used  in  this  report  lo  illustrate  the  various  model  inputs 
and  outputs  are  taken  from  the  analvsis  of  Lock  and  Dam  1  area.  The  tables 
and  alphameric  maps  employed  are  representations  of  the  modeled  area;  each 
symbol  or  figure  represents  a  value  for  a  grid  node  (which  represents  an  area 
0.3  bv  0.3  mi,  or  0.8  by  0.8  km). 

To  provide  for  continuitv  in  modeling  the  entire  navigation  reach,  the 
models  were  designed  to  include  an  area  of  overlap  on  the  adjacent  model. 

Adjacent  models  were  overlapped  a  minimum  distance  equivalent  to  <3  mi  (9.7 
km).  This  overlap  aided  in  the  identification  of  errors  associated  with  model 
boundarv  conditions  and  enabled  the  preparation  of  a  complete  suite  of  data 
for  tiie  navigation  reach.  As  the  models  for  adjacent  areas  were  annlvzed,  the 
data  developed  for  areas  common  to  each  model  were  examined  and  compared  to 
determine  the  extent  of  boundarv  effects.  Model  boundaries  parallel  to  the 
river  were  placed  at  a  distance  far  enough  from  the  river  so  that  the  ot Sects 
of  r iver- induced  water-level  changes  would  not  extend  to  the  boundaries. 

Nonstoadv  State 

Nonsteadv-state  analvscs  for  the  investigation  were  made  bv  using  three 
digital  programs  called  SFFKRMOCK,  DATE,  and  I1YDKOC  (Reed  and  others,  1976), 
which  were  developed  particularly  lor  this  studv.  SFI’ERMOCK  was  designed  to 
simulate  transient  stress  anil  response  in  a  ground-water  t  low  svstem  that 
includes  a  water  table  in  the  confining  laver  above  an  artesian  aquiter.  The 
model  incorporates  all  the  components  of  stress  in  the  t  low  field.  SiTERMtHlK 
mode  1 s  three  component  lavers:  a  so i l -mo i st ure- account i ng  component,  a  vertical- 
flow  component,  and  a  hor i zont a  1- f 1 ow  component  (t  ig.  3).  DATE  assigns  calendar 


Table  1. — Specifications  for  lock  and  dam  arrangements  studied  in  the 

investigation 
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Figure  5. --Relation  between  soil-moisture-accounting,  vertical-flow,  and 
horizontal-flow  components  of  SUPERMOCK  program.  (From  Reed 
and  others,  1976.] 
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dates  to  da  La  computed  at  observation  nodes  in  SI'  I’FRMi  K'K  and,  for  calibration, 
makes  comparisons  between  coi:i|)uted  potent  iometr  ic  and  water-table  level-,  1  ro- 
SI  PilkMlKlK  ami  observed  :  i  e  Id  values.  HYbROll  produces  livdrogranhs  ns  i  n.'.  out  sat 
1  rom  DAl'K.  liie  usr  oi  output  i  rovi  !1.\TF  and  HYDROO  is  disnissei;  later  in  t  h  i 
report  under  "iialilir.it  ion  ol  Moil,  Ls." 

For  ease  in  parameter  modi!  iealion  and  lor  adequate  model  in. ■  control, 
data  components  are  read  into  the  model  bv  discrete  subareas,  each  con  t  a i n  i  as 
one  observation  wr 1 1  at  which  control  data  bad  been  collected.  Is,  .am: isuru- 
t  ion  ol  these  subareas  was  determined  using  the  Til  iessen  polygon  a-tand.  '»  it  a 
that  arc'  entered  into,  or  computed  bv,  SlM’KR.MObK  in  this  manner  are:  the 
hydraulic  cominc t i v i t  ies  ot  the  upper  coniining  layer,  aquifer  storage  co.  :::- 
C  ients,  spec  i  f  i  c-v  i  e  1  d  values  ot  tlu1  confining  Laver,  and  evapo  t  ran  s  n  i  r.i  t  :  iitiy 
The  values  o!  these  parameters  were  determined  by  cal  i  brat  ion.  In  order  to 
maintain  control  over  calibration  cause  and  etlect,  t  he  value  designat'd  for 
each  ol  these  parameters  at  the  control  points  was  assigned  to  ill  nodes 
within  each  polvgon  (subarea)  in  the  model  area. 

l'iie  so  i  1 -mo  1st  ure-account  i  ng  component  in  Shi’FRMi'CK  is  a  parametric 
rainfall  accretion  mode  L  in  which  the  parameters  have  physical  signi:  icance. 
This  component  computes  changes  in  so  i  1-mo  i  s L u re  storage,  and  recharge  to  and 
discharge  from  the  zone  of  aeration  to  the  water  table.  Seven  parameter^  u-ed 
in  the  soil-mo isture-aocount ing  procedure  ho  1 p  define  tile  hydraulics  of  t  he 
soil  as  related  to  infiltration,  storage,  and  drainage.  file  values  o:  t  he  -v 
parameters  were  chosen  arbitrarily  by  a  t r i a  1-and-error  procedure  in  which 
infiltration  was  computed  based  upon  the  value  of  the  soils  parameters  and 
daily  precipitation  and  evaporation.  Plausibility  limits  for  the  parameters 
are  defined  in  the  soils  literature.  Within  tiiese  limits,  the  values  of  the 
soils  parameters  were  adjusted  until  a  combination  was  found  that  produced 
reasonable  infiltration  rates  for  the  types  of  soils  found  in  the  Red  River 
Valley.  These  seven  parameters  are: 

-  SMS  IN  - 

This  parameter  defines  the  initial  value  for  surface-mo i sturo  storage,  in 
inches.  Surface-moisture  storage  (SMS)  is  carried  bv  the  model  in  an 
array  containing  values  for  SMS  at  each  node  in  the  grid.  Tn  the  I irst 
Lime  step  of  tlu>  model,  each  member  of  this  arrav  is  sot  equal  to  SMS1N. 
For  the  Rod  River  models,  the  value  used  in  each  1 ock-and-dam  area  was 
1.0  in.  ( 2  > . 4  mm) . 


-  KSAT  - 

Tii  is  parameter  del  ines  tin-  saturated  hvdraulir  conductivity  lor  soil,  in 
inches  per  dav.  For  the  Red  River  models,  a  value  oi  10.0  in/d  ( J > s 
mm/d)  was  used.  This  value  was  within  plausible  limits  and  sicmed  to 
produce  the  best  results  based  upon  observed  dita. 


-  URN  - 

This  parameter  defines  the  maximum  drainage  rate  for  soil,  in  inches  per 
dav.  it  controls  the  amount  of  infiltration,  or  recharge,  to  the  water 
table  when  an  excess  in  soil  moisture  is  available.  A  value  of  ID.  0  in/d 
(254  mm/d)  was  used  in  the  Red  River  models. 


-  SWF  - 

This  parameter  defines  the  suction  (tension)  of  tiie  soil  at  field  capacity, 
in  inches.  The  value  used  in  the  Red  River  models  was  120  in.  (3,050 
mm).  This  is  a  typical  value  for  soils  in  the  project  area  and  was 
obtained  from  the  soils  literature. 


-  Rt'.F  - 

This  parameter  defines  the  ratio  of  wi  1 1 i ng-point  tension  to  tension  at 
field  capacity  (dimensionless).  The  value  used  in  the  Red  River  models 
was  40.0.  This  value  also  was  obtained  from  the  soils  literature  and  is 
a  typical  value  for  the  project  area. 


-  SMSM  - 


This  parameter  defines  the  maximum  amount  of  water,  in  inches,  that  can 
be  held  in  surface-moisture  storage.  The  value  for  SMSM  was  obtained  bv 
a  calibration  process  in  which  observed  hvdrographs  at  control  wells  were 
compared  with  computed  hvdrographs  at  the  same  locations.  A  value  of  1.0 
in.  (25.4  mm)  for  this  parameter  was  used  in  the  Red  River  models. 


-  XNORM  - 

This  dimensionless  parameter  defines  the  limits  of  the  recharge  rate.  it 
was  set  to  3  in  all  models  of  the  Red  River.  This  value  allows  the 
recharge  rate  to  range  from  zero,  for  SMS  y  ().5x(SMSM),  to  0.  1  5x(DR.V) , 
for  SMS=SMSM. 

The  value  of  each  of  these  parameters  was  held  constant  for  the  entire  model 
and  was  entered  to  SUPERMOCK  on  a  data  input  card. 

The  stress  on  the  soil-moisture-account ing  component  is  the  daily  differ¬ 
ence  between  precipitation  and  potential  evapotranspiration  which  is  input  to 
SUPERMOCK  on  cards.  When  the  stress  is  positive,  infiltration  to  soil  moisture 
is  computed  as  a  function  of  precipitation  in  excess  of  evapotranspiration, 
the  amount  of  moisture  already  in  storage,  and  the  hydraulic  properties  of  the 
soil.  Infiltration,  or  positive  downward  flux,  is  computed  bv  the  model, 
using  a  modified  version  of  a  routine  from  a  model  bv  Dawdv,  I.ichtv,  and 
Bergmann  (1972,  p.  B5-B8).  This  routine,  which  uses  4-minute  rainfall  periods, 
was  modified  to  correspond  to  the  1-day  rainfall  periods  used  in  this  model. 
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Overland  runoff,  or  infiltration  residual  computed  in  the  routine  of  Dawdv, 
Lichtv,  and  Bergmann  (1972),  was  dropped  from  the  so i 1-mois turn-account ing 
procedure  in  SUPERMOCK.  Due  to  the  1-dav  rainfall  period,  it  was  necessary  to 
impose  an  upper  limit  (SMSM),  as  previously  mentioned,  on  soil-moisture 
storage  because  redistribution  of  moisture  occurred  only  once  each  day.  The 
value  of  ttiis  limit  used  in  the  Red  River  models  was  l  in.  (25.4  mm).  Because 
the  surficial  material  of  the  Red  River  alluvium  is  generally  fine  grained,  a 
limit  of  soil-moisture  storage  of  1  in.  (25.4  mm)  is  reasonable.  Evapotrans- 
piration,  or  negative  stress,  is  subtracted  from  soil-moisture  storage  up  to 
the  amount  of  water  available.  When  soil  moisture  is  reduced  to  zero,  evapo- 
transpiration  is  derived  from  ground-water  storage  in  the  water-table  zone  in 
the  confining  bed  until  soil  moisture  is  replenished  from  rainfall. 

The  vertical-flow  component  in  SUPERMOCK  computes  the  elevation  of  the 
water  table  in  the  fine-grained  material  above  the  aquifer  a.-,  a  function  of 
the  elevation  of  the  water  table  in  the  preceding  time  step,  the  elevation  of 
the  potentiometric  surface,  and  recharge  from  the  soil-moisture  zone.  By  use 
of  this  water-table  elevation,  flow  to  or  from  the  aquifer  can  be  determined 
and  used  by  the  horizontal-flow  component.  SUPERMOCK  computes  the  redistribution 
of  soil  moisture  (recharge)  to  the  water  table  as  a  decaying  exponential 
function  of  soil  moisture  throughout  the  range  from  1  to  0.5  in.  (25.4  to  12.7 
mm).  For  soil  moisture  less  than  0.5  in.  (12.7  mm),  SUPERMOCK  sets  recharge 
to  the  water  table  to  zero.  Initially,  the  model  takes  evapctranspiration 
from  soil  moisture  and  then  from  ground-water  storage  in  the  upper  confining 
layer  after  soil  moisture  is  depleted.  The  limit  on  evapotranspiration  from 
ground  water  is  the  steady-state  rate  of  upward  movement  of  water,  as  deter¬ 
mined  by  the  method  of  Ripple,  Rubin,  and  van  Hylckama  (1972).  ATMOFLUX,  a 
peripheral  data-preparation  program  developed  for  the  investigation,  was  used 
to  compute  these  data.  ATMOFLUX  uses  a  method  requiring  a  specified  relation 
between  unsaturated  hydraulic  conductivity  and  soil  suction  (Ripple  and  others, 
1972,  p.  AA,  eq.  10).  Two  parameters  of  this  specification,  n,  an  integer 
soil  coefficient,  and  Sh,  soil  suction  at  which  the  unsaturated  conductivity 
is  one-half  the  saturated  conduct ivity,  are  used  to  express  the  limiting 
steady-state  evapo transp i  ra t  ion  in  a  nondimens ional  form.  Values  of  n_, 
ranging  from  2  for  clays  to  5  for  sands,  and  values  of  S^,  ranging  from  1  for 
sands  to  2  for  finer  materials,  were  used  in  this  study.  Output  from  ATMOFLUX 
includes  punched  cards  containing  values  of  evapctranspiration  divided  by 
saturated  hydraulic  conductivity  for  depths  to  the  water  table  ranging  from  1 
to  30  ft  (0.3  to  9.1  m)  for  four  ranges  in  hydraulic  conductivity  associated 
with  each  soil  coefficient,  n.  These  pinched  cards  are  used  as  input  to 
SUPERMOCK.  The  actual  limiting  rate  ol  evapctranspiration  used  by  SUPERMOCK 
was  obtained  by  multiplying  the  computed  upward  rate  associated  with  depth  to 
the  water  table  at  a  particular  time  bv  the  saturated  hydraulic  conductivity 
of  tite  upper  segment  (HCU)  of  the  upper  confining  layer  in  a  particular  subarea. 
The  method  of  Ripple,  Rubin,  and  van  Hylckama  (1972)  assumed  hare  scil  anil 
moisture  transport  to  the  land  surface.  Practical  1 v  all  the  Red  River  project 
area  is  covered  by  vegetation.  Therefore,  moisture  transport  was  calculated 
to  the  base  of  the  root  zone. 
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The  hor  izontal-f  low  component  in  SUPF.RMOCK  computes  the  transient  eleva¬ 
tion  of  the  potent iometric  surface  in  the  aquifer.  In  the  Red  River  models, 
the  stresses  on  the  aquifer  that  were  simulated  included  the  imposition  of 
time-variant  stream  stages  for  tiie  main  stem  of  the  Red  River  and  its  major 
tributaries  and  accretion,  which  is  computed  by  SUPERMOCK  as  a  function  of  the 
water-table  elevation.  Where  a  computed  water  table  does  not  exist,  the  model 
uses  infiltration,  or  recharge,  from  the  so i 1 -moisture  zone  as  accretion  to 
the  aquifer. 

The  time-step  increment  used  in  the  nonsteady-state  analyses  of  the  Red 
River  models  was  10  days.  Time-variant  stream-stage  and  climatic  data  were 
used  as  input,  and  the  potentiometric  surface  and  water-table  elevations  at 
each  node  in  the  grid  were  computed  for  each  time  step. 

Calibration  of  the  nonsteady-state  model  was  based  upon  preconstruction 
stream  stages  and  comparisons  of  computed  and  observed  hydrographs  at  observa¬ 
tion  wells.  After  calibration,  the  model  was  used  to  compute  postconstruction 
elevations  of  the  potentiometric  surface  and  water  table.  Postconstruction 
output  was  based  upon  the  imposition  of  postconstruction  stream  stages  on  the 
main  stem  of  the  Red  River.  The  availability  of  the  time-varying  elevation  of 
the  water  table  allowed  the  computation  of  average  depths  to  the  water  table  for 
specific  periods  of  interest  requested  by  the  Soil  Conservation  Service. 


Steady  State 

Steady-state  projections  of  the  postconstruction  potentiometric  surface 
in  the  Red  River  alluvial  aquifer  were  made  using  techniques  developed  during 
similar  studies  in  the  Arkansas  River  valley  (Bedinger  and  others,  1970). 
During  the  Arkansas  River  study,  these  techniques  were  applied  to  analog 
modeling.  For  the  Red  River  investigation,  these  techniques  were  incorporated 
into  a  digital  model  called  GWFLOW  (Bedinger  and  others,  1973).  GWFLOW  is  a 
two-dimensional  representation  of  an  aquifer. 

The  principal  data  needs  of  the  GWFLOW  model  for  use  in  steady-state 
analysis  are  transmissivity  of  the  aquifer,  the  ratio  of  change  in  evapo- 
transp irat ion  to  change  in  aquifer  head  (AFT/AH),  change  in  stream  stages,  and 
thickness  and  hydraulic  conductivity  of  streambed  material.  To  determine  the 
change  in  head  at  any  point  in  the  aquifer  resulting  from  a  change  in  river 
stage,  the  initial  potentiometric  surface  on  the  stream  boundaries  is  the 
change  in  river  stage  and  is  zero  at  all  other  nodes  in  the  aquifer. 

In  the  steady-state  models  of  the  Red  River  alluvial  aquifer,  trans¬ 
missivity  was  varied  over  the  modeled  area,  and  AF.T/AH  was  entered  as  varying 
by  discrete  subareas.  The  method  used  to  determine  values  of  AET/AH  is  dis¬ 
cussed  later  under  "Preparation  of  Digital-Model  Input"  and  "Calibration  of 
Models."  Stress  on  the  models  was  imposed  at  appropriate  stream  nodes  as 
changes  in  stream  stage  from  preconstruction  to  pos tcons t rue t ion  conditions. 
The  direct  effects  of  changes  in  stage  for  streams  with  partial  hydraulic 
connection  were  simulated  by  applying  nonuniform  streambed  thickness  and 
holding  the  hydraulic  conductivity  of  the  streambed  material  constant.  The 
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values  of  AET/AH  have  a  definite  controlling  effect  on  the  magnitude  of  change 
in  the  po ten tiome trie  surface  and  on  the  area  of  influence  of  stream-stage 
change . 

Time-step  increments  for  GWFLOW  were  based  on  computation  times  entered 
on  cards.  The  computation  times  used  in  the  Red  River  models,  which  were 
those  that  were  recommended  for  GWFLOW,  ranged  from  0.00110  to  40,000  days  in 
logarithmic  increments.  Although  analyses  indicated  that  most  of  the  water- 
level  changes  had  taken  place  in  the  first  2-1  years,  computation  Limes  were 
extended  to  40,000  days  to  insure  complete  equilibrium.  Primary  output  from 
the  models  consisted  of  changes  in  the  potent iometr  ic  surface  at  each  node  in 
the  9.5-mile  (0.8-km)  grid.  This  output  was  used  to  contour  changes  in  the 
potent iometr ic  surface  in  the  aquifer  resulting  from  an  increase  in  river 
stage . 


PREPARATION  OF  DIGITAL- MODEL  INPUT 

Preparation  of  input  data  for  use  in  the  GWFLOW  and  SUPERMOCK  models 
involved  the  collection  and  manipulation  of  field  data.  Some  of  the  data 
required,  and  also  the  data  format,  are  common  to  both  GWFLOW  and  SUPERMOCK. 
However,  because  of  the  greater  complexity  of  the  SUPERMOCK  model,  more  detailed 
and  varied  types  of  input  were  required  for  it  than  for  the  GWFLOW  model. 

Several  data-preparation  computer  programs,  hereinafter  termed  "peripheral 
programs,"  were  developed  during  the  investigation  to  process  data  required  by 
the  models.  These  programs  will  be  discussed  in  the  following  sections. 

Source  listings  and  data-input  requirements  of  these  peripheral  programs  are 
included  as  attachments  at  the  end  of  this  report. 

Some  of  the  data  read  into  GWFLOW  were  dependent  upon  parameter  values 
determined  during  the  calibration  of  the  nonsteady-state  model.  Therefore, 
nonsteady-state  analyses  for  each  lock-and-dam  area  were  made  before  the 
corresponding  steady-state  analyses  for  that  area.  For  purposes  of  discussion, 
preparation  of  data  for  the  two  models  will  also  be  discussed  in  that  order. 


Nonsteady-Sta Le  Model 

Varied  types  of  data  were  prepared  for  entry  into  the  nonsteady-state 
model  in  order  to  adequately  define  the  flow  field.  Most  of  this  input  is  in 
the  form  of  alphameric  maps  that  are  representations  of  the  modeled  area. 

Many  of  these  maps  are  outputs  from  the  peripheral  programs  mentioned  previously. 
The  primary  data  input  to  the  model  are  depicted  in  the  generalized  flow  chart 
in  figure  6. 

Root  Depth 

Root  depths  of  vegetation  are  key  factors  required  by  SUPERMOCK  in 
determining  the  effective  depth  to  the  water  table  for  computation  of  evapotrans- 
piration.  Evapotranspiration  is  modeled  as  depleting  the  moisture  content  in 
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the  soil  Liver  between  land  surface  and  the  base  ol  the  root  zone.  Upward 
t low  from  tile  water  table  occurs  as  a  response  to  this  surfirial  depletion. 

Trapping  patterns  to  the  nearest  40-acre  ( 1 b2 , 00()-m~ )  plot  and  the  effec¬ 
tive  root  depths  tor  tlie  various  types  of  vegetal  cover  were  determined  bv  the 
Soil  Lonservat  ion  Service.  This  information  was  based  upon  1971  cropping 
patterns  that  were  assumed  to  be  representative  of  the  project  area  for  the 
calibration  period. 

The  root-depth  data  tor  each  lock-and-dam  area  were  entered  into  SUPERMOCK 
in  the  torm  ot  an  alphameric  map  on  cards.  The  various  tvpes  of  vegetal 
cover,  their  associated  root  depths,  and  the  svmbols  representing  those  depths 
are  tabulated  below: 


Vegetal  cover 


Root  depth 
(feet) 


Map  symbol 


Cotton -  2.3  C 

Soybeans -  2.3  S 

Pasture -  2.5  P 

Orchards -  5.0  0 

Woodlands -  5.0  W 

Uplands -  3.0  U 

Urban  areas -  2.5  E 


An  example  of  an  alphameric  root-depth  map  is  shown  in  figure  7.  (All 
examples  are  for  Lock  and  Dam  3  area.) 


Land-Surface  Elevation 

Fhe  elevation  of  the  land  surface,  in  feet  above  Mean  Sea  Level  Datum  of 
1929  (now  referred  to  as  National  Geodetic  Vertical  Datum  of  1929),  was  used 
in  the  nonsteady-state  models  as  a  reference  point  for  determining  (1)  the 
depth  to  the  water  table,  (2)  the  relation  of  the  potentiometric  surface  to 
land  surface,  and  (3)  the  elevation  of  the  top  of  the  aquifer. 

Land-surface  elevations  were  obtained  from  two  sources:  instrument  levels 
and  topographic  maps.  The  land-surface  elevation  at  each  of  the  observation 
wells  was  determined  by  instrument  and  assigned  to  the  node  nearest  the  well. 

At  all  other  nodes  in  each  of  the  lock-and-dam-area  models,  these  data  were 
picked  from  topographic  maps.  The  appropriate  set  of  data  was  read  into 
SUPF.RMOCK  for  each  lock-and-dam  area  in  the  form  of  a  numeric  map.  Land- 
surface  elevations  at  each  node  were  estimated  to  the  nearest  foot. 

Topographic  map  coverage,  including  map  contour  interval,  is  shown  in 
figure  8. 
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Initial  Po  t  on  t  iome  I  r  i  e  Surface  and  Water  Table 


Tlu'  elevation  of  thi'  water  table  .mil  pot  ent  i  omol  r  i  e  niri  .uv  i  s  reqii  i  r<  «i 
bv  Sl'PKRMOCK  .is  tlu'  starting  point  tor  computations.  In  tin  t  irst  t  i  :u.  st.-p 
in  tin'  moili'l,  tin-  elevation  ot  the  water  table  is  set  eipial  to  tin-  elevation 
ot  tin'  pot  i'll  t  i  omet  r  i  e  surl  are  at  eor  respond  i  ng  nodes.  liieretoro,  the  data 
input  that  must  hi'  prepared  is  the  initial  potent  iiinietrii'  surtaee. 

In  tin'  Red  Kiver  models,  the  average  pot en t i omet r i v  surtaee  tor  the 
period  ot  record  was  used  to  represent  the  init  ial  p recoils t rue t i on  potent  i o- 
tlietric  surtaee.  The  elevation  ot  the  average  potent  iomet  lie  surtaee  in  the 
aquifer  was  based  on  a  minimum  ot  4  vears  ot  record.  these  data  were  collected 
I  ruin  the  joint  Geological  Survev-Soil  Conservation  Service  observation-well 
network  in  tlu*  valley.  Measurements  were  made  monthly  in  the  iV)-well  network. 
Water-level  measurements  at  each  observat ion  well  in  a  part  icuiar  lovk-and-.lam 
area  were  averaged  on  a  time-weighted  basis  using  a  digital  urogram  called 
AVERAGE,  which  was  developed  tor  this  purpose.  l'he  onlv  data  required  hv  tin 
AVERAGE  program  are  water  levels  and  corresponding  dates  ot  measurement  at 
each  observation  well.  A  program-source  listing,  containing  input  requirements 
and  formats,  and  an  example  ot  program  output  are  included  in  attachment  A. 
l’he  average  values  determined  t  rotn  this  procedure  were  plotted  and  manually 
contoured  to  obtain  t  lie  elevation  ot  the  p  recons  t  rue  t  i  on  pot  oil  t  iomet  r  i  e 
surface  in  that  lock-and-dam  area.  I'lie  resulting  average  potent  iomet  r  ie 
surface  represents  a  hypothetical  dvnam i c-equ i 1 i hr i urn  condition  of  luad  in  the 
aquifer  for  preconst ruct ion  conditions.  thorn  the  contour  map,  the  elevation 
of  the  potent iomet r ic  surface  was  picked  for  each  node  in  the  grid  covering  a 
lock-and-dam  area.  These  values  were  coded  into  a  numeric  map  containing 
elevations  to  the  nearest  toot  at  each  node.  I'lie  map  was  converted  into  data 
cards  that  were  used  as  input  to  SITEKMOCK. 
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Observed  1’otent  iomet  r  ic  Surtaee  and  Water-table  Klevat  ions 

For  purposes  of  calibration,  observed  levels  ot  tin*  water  table  and  the 
potent iomet r ic  surface  were  compared  with  corresponding,  values  computed  hv 
SUPERMOCK.  Tile  comparisons  were  made  in  the  DATE  program  (Reed  and  otlur-., 

19  76)  that  was  run  in  sequence  with  the  StTKRMOCK  model.  l'he  observed  data 
used  in  DATE  consisted  of  the  spring  "high"  and  tall  "low"  water  table  and 
po ten t  iome t r i c  levels  tor  one  or  more  vears. 

The  observed  potent  iomet  r  ic  levels  for  the  high  in  the  spring,  and  the  low 
in  the  fall  of  specified  vears  were  read  into  DATE  as  exact  valuer.  However, 
because  the  position  ot  the  water  table  at  some  sites  known  onlv  within  a 
certain  range,  several  input  tormat  options  are  allowed.  Water- t able  value' 
mav  be  entered  as  being  greater  than  or  less  than  a  given  value,  a  .  being 
within  a  closed  range,  as  an  exact  value,  or  as  being  unknown. 

An  example  of  a  calibration  table  produced  bv  DA  IT  and  a  disco. -ion  ot 
the  use  of  the  data  are  given  in  the  section  "Ca 1  i  h ra t  i on  and  V>  ■  r  i  t  i  cut  i  on  ot 
the  Nonsteady-State  Model."  Observed  data  are  printed  in  the  table  aoeorJing. 
to  the  format  in  which  tliev  were  entered  ;o  DA  IT.  Use  ot  the  observed  data 
for  comparisons  with  computed  data  was  invaluable  in  the  •alibr.it  ion  imocoss. 
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transmissivity  of  the  alluvial  aquifer  in  the  study  area  ranges  from 
11,000  ft~/d  (279  to  1,290  m-/d)-  These  values  were  determined  at 
sites  hv  analysis  of  pumping-test  data  and  hv  analysis  of  aquifer 
to  river-stage  fluctuations.  These  data  were  extrapolated  to  other 
the  valley  by  developing  relationships  between  hydraulic  conductivity 
cle  size  at  the  pumping-test  sites  and  extending  these  values,  on  t  In¬ 
grain-size  relationships  and  thickness,  to  test-hole  sites. 


Pumping  tests  conducted  by  the  Geological  Survey  as  part  of  earlier 
studies  ol  the  alluvium  (Newcome,  19b0)  provided  values  of  transmissivity  at 
six  locations  in  the  valley.  Transmissivity  values,  determined  from  these 
tests,  ranged  from  1,300  to  13,000  ft^/d  (492  to  1,210  m"/d).  The  hydraulic 
conductivity  ranged  from  1)0  to  160  ft/d  (40  to  49  m/d). 


Approximately  110  samples  of  aquifer  material  were  collected  from  test 
holes  and  analyzed  for  hydraulic  conductivity  and  particle  size  during  the 
investigation.  From  these  analyses,  a  relationship  was  developed  between 
hydraulic  conductivity  and  particle  size,  using  Lin-  method  ol  Johnson  and 
bedinger  (19b7).  From  this  relationship,  an  average  value  ot  hydraulic 
conductivity  was  developed  for  the  alluvial  aquifer.  Genuine' t i v  i  t v  values 
obtained  hv  this  method  were  compared  with  those  determined  I rom  pumping 
tests.  From  these-  analyses,  an  average  value-  eif  hvdraulic  con  due  t  i  v  i  t  v  e>f  147 
ft/d  (4i  m/d)  was  determined  lor  the  alluvial  aquifer.  This  value-  was  clu-i-ke-d 
at  several  lex-ations  near  the  rive-r  hv  using  the-  R I VKR- 1  NTH'GKD  I'l.l'CTl’AT  I ONS 
cemipute-r  program  (Be-eiinge-r  and  others,  197  1).  The-  transmissivity  at  e-ach  ol 
the-  te-st— hole'  sites  was  then  computed  hv  multiplying  the-  ave-rage-  oetneluc  t  i  v  i  t  \ 
hv  t  in-  thickne-ss  ol  aquife-r  mate-rial  note*d  in  the-  test -hole  logs. 

1  ran  sm  i  s  s  i  v  i  t  v  va  1  ues  for  t  In-  terrace  tie- posits  were-  estimates!  using 
t  h  i  e'kin-sse-s  obtained  fremi  logs  of  te-st  holes  in  the'  eie-posits.  Tiu-  average' 

,'ivel  ran  !  ie-  coneluc  t  i  v  i  t  v  was  assumed  t  e'  be  147  ft/d  (4>  m/el).  I'errae-e  el.-posits 
were'  assigne-d  transmissivity  value's  where  t  lie'v  are  areal  1  v  extensive'  aiul  are' 
e'eilis  i  dere-el  t  e'  he  h  vel  rail  1  i  e-a  1  1  V  e-einne'C  t  e'd  with  tile'  alluvial  aquiteT. 

1'iie  I  orina  t  i  imis  eH  Te'rtiarv  age',  vliie-li  underlie'  the-  alluvium  and  le'rm  the 
upl.nnls  borde- r  i  tig  the'  valley,  are'  e-e'mposenl  primarily  i>!  silt  ami  e' lav  ami  an- 
relative!',  i  inperine.ih  1  <•  compare-el  with  the  alluvial  aquile-r.  l.-.t  imat  esi  t  i  aib- 
m  i  -;s  i  i  t  i  e's  i  ei  r  sami  units  in  t  hc-so  leirmalieins  range-el  t  re«m  >t‘  t  e>  /  I  t  •  A! 

(4./  tei  ni  nr  /ei)  in  areas  where-  thee  are-  in  hvelraulie  cemiie-e'l  imi  with  tin' 
alluvia!  aqu  iter.  These  estimates,  wore-  based  upeui  geo  1  og  i  e-  .me!  pump  iuv-t  es.t 
eiala  e  ei  1  1  in' t e  ei  during  esarlier  slinlies  (Se'wceimi-,  ll,n0i. 


Alter  transmissivity  value",  h.iel  lieu-n  nlotteel  ami  evn  t  mi  i  s  J  i  t>r  tin 
pro'll'.  I  ate. i,  alphameric  maps  were  pt  epar<-el  :  or  <  ach  link-ami-. lam  area;  am! 
tin  dat  i  wer.  piimhe.l  on  eards  loi  input  l  .>  I  in-  -uni.  1  .  \n  .-va'-ip  I  •  .»!  an 
a  1  chanter  i  c  t  I  all  -.m  i  ■■  s  j  v  j  t  \  map  I  rom  tin  I  udv  ami  e  xp  I  .mat  ion  o!  -.vmluil  are 
.  howil  in  I  i  ..pi  r.  '* . 
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Figure  9. --Example  of  alphameric  transmissivity  map 


Additional  chocks  were  made  on  the  modeled  transmissivity  values  during 
calibration  of  the  nonsteadv-state  models.  However,  onlv  minor  adjustments 
were  made,  and  the  maps  were  used  virtually  as  initiallv  prepared  in  both  the 
steadv-  and  nonsteady-state  models. 

Conductivity  of  the  I'pper  Confining  Laver 

Movement  of  water  to  or  from  the  alluvial  aquifer  takes  place  through  tin- 
upper  cent  ining  layer,  which  overlies  the  aquifer  nearly  everywhere  in  the 
valley.  The  upper  confining  laver,  which  is  composed  of  a  heterogeneous 
sequence  of  clav,  silt,  and  sand,  ranges  in  thickness  from  a  few  feet  to  140 
tt  (41  m).  Movement  ol  water  through  the  upper  confining  laver  was  modeled  as 
being  one-dimensional  vertical  1  low.  To  provide  for  greater  flexibility  in 
mod  ling  the  vertical-flow  component,  the  upper  confining  layer  was  modeled  as 
two  segments;  one  segment  extending  from  the-  base-  of  the  root  zone  to  tin- 
water  table  and  the  other  extending  from  the  water  table  to  the  top  of  the 
aqu iter. 

Both  the  upper  and  lower  segments  were  assigned  values  of  hydraulic 
conduct  iv  i  tv ,  designated  HCll  and  11CL,  respectively.  These  values  were  entered 
in  SITKRMOCK  bv  discrete  subareas— each  subarea  having  a  unique  value  for  HCL 
and  He'. L .  Kacli  node  within  a  subarea  was  assigned  the  same  value  for  HCl'  and 
HCL.  Initiallv,  1ICL  and  HCL  values  in  a  particular  subarea  were  set  equal  to 
the  same  value.  This  value  represented  the  harmonic  mean  of  the  conduc t i v i t i e 
for  materials  in  the  upper  confining  layer  in  that  subarea.  These  harmonic- 
mean  conductivities  were  computed  using  a  digital  program,  ATMOl'LUX,  shown  in 
attachment  B.  The  ATMOI-'Ll'X  program  uses  as  input  the  thickness  and  lithologic 
class  for  materials  in  the  upper  confining  laver.  Lithologic  data  lor  the 
upper  confining  laver  were  obtained  from  test-hole  logs.  The  scheme  used  in 
the  studv  for  associating  lithologic  class  and  hvdraulic  conductivity  is  shown 
in  figure  10. 

llvdraul  ic-conduc  t  i  v  i  t  v  values  ranging,  I  rom  1.0  to  1.0x10  ’  It/d  (O.o  to 
3xl0_fl  m/d)  were  selected  as  being  the  physical  plausibility  limits  within 
which  adjustments  could  be  made  to  the  vertical  hvdraulic  londuetivitv  ol  the 
upper  coni  ining  laver.  I'll  i  s  range  represents  the  conduc  t  i  ,>  i  t  v  ol  materials 
ranging  from  I ine  sand  to  dense  clav.  Because  ol  the  lateral  variability  of 
upper  alluvial  materials,  the  initial  conductivity  values,  as  determined  1 rom 
test-hole  logs,  are  not  necessarily  representative  ol  the  entire  area  as 
modeled.  Therefore,  the  onlv  constraints  on  adjusting  vertical  hvdraulic 
conductivity  values  during  calibration  was  to  remain  within  the  physical 
plausibility  1  im its. 

An  example  of  an  alphameric  map  and  the  accompanying  table  de' ining  the 
value  of  IICU  and  HCL  for  each  subarea  of  a  lock-and-dam  area  are  shown  in 
t  igure  1 l . 
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-  Example  of  alphameric  map  of  vertical  hydraulic  conductivity  and 
explanation  of  symbols 


Relation  of  Evapotransp i rat  ion  to  Depth  to  Water 


Because  the  entire  valley  is  covered  hv  vegetation,  the  removal  of  water 
by  evapot ranspi rat  ion  is  not  at  the  land  surface  but  is  at  the  base  of  the 
root  zone  in  the  fine-grained  layer.  To  determine  the  rate  of  evapotrans- 
piration  from  the  root  zone  for  different  depths  to  water,  a  function  expressing 
the  relationship  between  dimensionless  evapotransp irat ion  and  depth  to  water 
below  the  root  zone  (GWETO)  is  used  bv  the  model.  Values  of  the  GWETO  function 
were  computed  by  the  ATMOFLl’X  program  (attachment  F)  and  were  entered  to  the 
model  on  cards.  GWETO  includes  four  different  functional  relations  between 
evapot ranspi rat  Ion  and  saturated  hydraulic  conductivity.  The  values  of  the 
GWETO  function  for  the  four  ranges  in  hydraulic  conductivity  and  for  depths  of 
from  1  to  30  ft  (0.3  to  9.1  m)  to  the  water  table  are  shown  in  figure  12.  The 
appropriate  relation  is  chosen  during  program  execution  based  on  the  value  of 
HCU.  Tiie  value  of  evapotranspiration  is  computed  in  the  program  as  the  product 
of  GWETO  at  a  particular  depth  to  water  and  the  upper  hydraulic  conductivity 
(HCU).  A  detailed  discussion  of  the  determination  of  the  GWETO  function  is 
given  in  Reed,  Bedinger,  and  Terry  (1976,  p.  52). 


Figure  12.— Example  of  GWETO  functions  for  computation  of  evapotransoiration 


Thickness  of  Streambed  Material 


The  Red  River  and  its  tributaries  do  not  fully  penetrate  the  alluvial 
aquifer  at  all  places  along  their  channel.  The  fine-grained  material  that 
exists  beneath  the  stream  channels  in  places  retards  the  movement  of  water  to 
or  from  the  aquifer.  As  a  result,  for  preconstruction  conditions,  water 
levels  in  observation  wells  as  close  as  200  ft  (61  m)  to  the  streams  nay 
differ  bv  as  much  as  1  to  5  ft  (0.9  to  1.5  m)  from  stream  levels  during 
transient  conditions. 

In  SUPERMOCK,  grid  nodes  assigned  to  the  main  stem  of  the  Red  River  mav 
optionally  be  specified  as  fully  or  partially  penetrating  the  aquifer.  All 
tributary  stream  entries  are  assumed  to  be  partially  penetrating.  The  model 
requires  that  all  partially  penetrating  stream  nodes  be  assigned  a  streambed 
thickness . 

The  thickness  of  material  beneatli  the  stream  channels  was  not  known 
initially  except  through  qualitative  estimates  based  on  logs  of  test  holes 
near  the  stream  channels.  Therefore,  the  effective  thickness  was  determined 
from  analysis  of  SUPERMOCK's  response  to  different  thicknesses  as  indicated  by 
the  differences  in  the  computed  and  observed  potentiometric  surface  at  control 
wells  near  a  stream.  The  reasonableness  of  the  annual  accretion  to  the 
aquifer  necessary  to  maintain  a  computed  potentiometric  level  equal  to  the 
observed  level  at  those  wells  was  also  considered.  An  arbitrary  value  of 
thickness  was  assigned  to  each  node  in  the  model  that  represents  a  point  on  a 
stream  channel.  Maps  showing  streambed  thickness  were  then  prepared  for  each 
of  the  modeled  areas.  Separate  symbols  were  used  for  each  stream,  and  an 
arbitrary  value  of  thickness  was  given  to  each  symbol.  During  calibration, 
additional  symbols  were  introduced  where  needed  to  represent  different  thick¬ 
nesses.  'There  changes  were  not  required,  the  symbols  used  initially  were 
retained  for  ease  in  identifying  various  modeled  stream  channels.  For  reaches 
of  a  stream  where  zero  thickness  seemed  to  be  indicated  by  model  response,  a 
very  small  nonzero  value  was  assigned.  The  program  logic  in  SUPERMOCK  computes 
no  flow  through  the  streambed  if  a  zero  thickness  is  coded  for  the  node  in  the 
streambed-th ickness  map.  An  example  of  a  streambed-thickness  map  and  its 
accompanying  legend  are  shown  in  figure  13.  The  thickness  value  associated 
with  the  symbols  H  and  C  is  printed  as  zero  because  of  the  print  format  in 
SUPERMOCK.  The  value  is  actually  a  small  nonzero  fraction.  The  blank  in  the 
e\planation  indicates  a  nonstream  node  and  therefore  has  no  streambed  thick¬ 
ness  associated  witli  it.  The  3's  around  the  edge  of  the  model  indicate  a  no¬ 
flow  boundary. 

The  thicknesses  shown  on  the  maps  do  not  necessarily  indicate  the  physical 
thickness  of  fine-grained  material  at  a  given  location.  A  single  value  of 
5xl0"3  ft/d  ( 1 . 5x10“ 3  m/d)  was  used  in  the  model  as  the  hydraulic  conductivity 
of  the  fine-grained  material.  Therefore,  the  thickness  was  adjusted  to  obtain 
the  correct  ratio  of  hydraulic  conductivity  to  thickness  for  calibration. 

Also,  because  of  the  0.5-mile  (0.8-km)  grid  spacing  used  in  the  models,  any 
modeled  watercourse  is  effectively  0.5  mi  (0.2  km)  wide.  Thus,  the  modeled 
thicknesses  must  represent  the  flow  characteristics  through  streambed  materials 
in  generally  much  narrower  streams.  A  near-zero  thickness  of  streambed 
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Map  of  thickness  of  streambed  and  lakebed  material 

SYMBOL  TMlCWESS  (ft ) 
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Figure  13.-Example  of  alphameric  streambed-thickness  map  tor 
nonsteady-state  analysis. 
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material  is  an  indication  that  at  that  point  the  river  and  aquifer  are  in 
perfect  hydraulic  connection. 

Changes  in  the  preconstruction  streambed-thickness  map  required  in  the 
postconstruction  analysis  included  only  the  addition  of  nodes  reflecting 
postconstruction  changes  in  the  position  of  the  navigation  channel.  The  added 
nodes  were  assigned  the  same  thickness  value  as  adjoining  river  nodes. 


Specific  Yield  and  Storage  Coefficient 

1’lie  introduction  of  specific-yield  values  and  aquifer-storage  coefficients 
into  the  model  was  done  in  alphameric  form  bv  discrete  subareas  identical  with 
those  used  for  the  entry  of  HCU  and  HCL.  An  example  of  an  alphameric  map  of 
specif  ic -yield  values  and  storage  coefficients  and  the  explanation  for  each 
are  shown  in  figure  14.  The  scheme  for  applying  calibration  values  to  identical 
subareas  was  used  so  that,  during  the  calibration  process,  modifications  could 
be  made  to  the  values  represented  by  symbols  in  anv  one  subarea  without  substan¬ 
tially  affecting  adjacent  areas. 

Specific-yield  values  were  limited  to  a  plausibility  from  lxlO~“  to 
2xl0~l,  and  the  storage  coefficient  was  allowed  to  vary  from  1x10”^  to  1x10"  \ 
The  final  specific-yield  values  and  aquifer  storage-eoef f ic ient  values  wore 
adjusted  in  the  calibration  procedure  bv  a  trial-and-error  process  within 
these  limits. 


Precipitation  and  Potential  Evapotranspi rat  ion 

Dailv  precipitation  and  evapotransp  i rat  ion  data  are  required  bv  Sl'PKRMOCK. 
Climatic  information  used  in  preparing  these  data  was  obtained  from  National 
Weather  Service  stations  in,  or  near,  each  lock-and-dam  area.  Data  from 
Weather  Service  stations  at  Alexandria,  Natchitoches,  Westdale,  and  Shreveport 
were  used  in  Lock  and  Dam  2-5  areas,  respectively  (fig.  3).  Data  for  Lock 
and  Dam  1  area  were  taken  from  the  Jonesville  station,  which  is  about  20 
mi  (32  km)  north  of  that  area. 

Daily  precipitation,  in  inches,  was  taken  directly  from  Weather  Service 
records  and  coded  for  card  input  to  SUPERMOCK.  The  model  assumes  uniform 
distribution  of  precipitation  throughout  the  grid  area.  Therefore,  no  nodal 
specifications  were  required.  SUPERMOCK  required  that  the  precipitation  data 
begin  on  or  before  the  first  day  of  the  simulation  run  and  continue  through 
the  duration  of  the  period  analyzed. 

Daily  potential-evapotranspirat ion  data  were  not  directlv  available. 
Therefore,  a  computation  scheme  was  required  to  derive  the  data.  Potential 
evapot ransp  irat ion  is  the  combination  of  evaporation  from  the  ground  stir  taco 
and  transpiration  from  plants  when  there  is  complete  vegetal  coverage  and  soil 
moisture  is  adequate.  Potential  evapotranspiration  was  computed  bv  the  method 
of  Thornthwaite  (1948).  This  computation  scheme  was  incorporated  into  a 
digital-computer  program  called  POTEET,  which  was  modified  from  a  program 
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developed  by  E.  P.  Weeks  (written  commun. ,  1973).  The  principal  data  require¬ 
ments  of  this  program  are  minimum  and  maximum  daily  air  temperatures,  monthly 
average  temperatures  during  the  period  for  which  potential  evapotranspi rat  ion 
is  to  be  computed,  and  latitude.  A  source  program  listing  and  complete  data 
requirements  for  POTEET  are  included  in  attachment  C.  Primary  output  f rom 
POTEET  consists  of  punched  computer  cards  that  are  in  a  format  compatible  with 
input  requirements  for  SUPERMOCK. 


River  Stage 

Two  complete  sets  of  time-variant  stream-stage  data  for  the  Red  River  and 
its  major  tributaries  were  required  for  nonsteadv-state  analysis.  Preconstruc¬ 
tion  conditions  in  each  lock-and-dam  area  were  simulated  and  the  nonsteadv- 
state  model  was  calibrated  to  reproduce  observed  water-table  levels  and 
potentiometr ic-surf ace  elevations  at  control  wells.  After  successful  calibra¬ 
tion,  the  preconstruction  stages  were  replaced  by  time-variant  postconstruction 
stages,  and  production  runs  were  made  simulating  postconstruction  conditions 
in  the  flow  field.  Datum  for  all  stream-stage  data  used  in  the  nonsteadv- 
state  model  was  Mean  Sea  Level  Datum  of  1929. 

The  Corps  of  Engineers  provided  time-variant  preconstruction  and  post¬ 
construction  stages  on  the  main  stem  of  the  Red  River  for  the  period  December 
1967  to  September  1973.  These  data  consisted  of  sets  of  5-dav-average  stages 
at  approximately  2-mile  (3.2-km)  intervals  for  the  entire  reach  of  the  Red 
River  in  the  project  area.  Each  set  of  associated  stage  and  river-mile  data 
was  identified  by  a  sequence  number,  increased  by  5  for  each  set,  to  correspond 
to  the  time  (day)  on  which  the  average  stages  were  based.  The  preconstruction 
and  postconstruction  stages  comprised  two  separate  data  sets,  each  residing  on 
a  separate  7-track  magnetic  tape.  These  data  sets  were  transferred  to  9-track 
tapes  and  used  as  master  input-data  sets  for  the  creation  of  separate  lock- 
and-dam-area  main-stem  river-stage-data  sets,  as  needed. 

The  individual  lock-and-dam-area  sets  for  the  main  stem  of  the  Red  River 
were  created  by  use  of  a  digital  program  called  RIVCHANC.E,  developed  specifi¬ 
cally  for  that  purpose.  The  source-program  listing  of  RIVCHANGE  and  data- 
input  requirements  and  formats  are  included  in  attachment  D.  Input  requirements 
for  RIVCHANGE  include  the  following:  (1)  beginning  and  ending  sequence  numbers 
corresponding  to  the  beginning  and  ending  dates  of  a  period  of  time  encompassing 
the  period  to  be  analyzed  for  a  particular  lock-and-dam  area;  (2)  a  number 
equal  to  an  interpolated  sequence  number  within  the  period  specified  in  (1)  at 
which  computation  of  10-day  averages  is  to  begin;  (3)  the  length  of  time,  in 
days,  for  which  computation  of  10-day  averages  is  to  continue;  (4)  the  beginning 
and  ending  river  miles  in  a  particular  lock-and-dam  area;  and  (5)  grid  nodes 
and  associated  river  miles  at  which  10-day-average  river  stages  were  desired. 
Node  designation  and  associated  river  mile  were  determined  manually,  beginning 
at  the  downstream  end  of  the  model  and  proceeding  upsteam  sequentially  to  the 
upstream  end  of  the  model  area. 

RIVCHANGE  was  designed  to  interpolate  in  time  and  space  and  compute  10- 
day-average  stages  at  specified  river  miles  associated  with  river-stage  nodes 
in  the  model  of  a  particular  lock-and-dam  area.  The  program  first  located  the 
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specified  time  period  within  the  master  data  set  and  determined  the  reach  of 
the  river  to  be  analyzed.  The  spatial  interpolation  was  based  on  river  miles 
and  the  temporal  interpolation  was  based  on  sequence  numbers  and  associated 
calendar  dates.  As  enough  daily  data  became  available  from  the  interpolation, 
RIVCHANGE  began  computing  10-day-average  stages,  beginning  with  the  dav 
designated  bv  the  beginning  sequence  number  for  computations,  and  continuing 
for  tiie  number  of  days  specified. 

Output  from  RIVCHANGE  consisted  of  10-dav-average  river  stages  associated 
witli  specified  grid  nodes.  Each  set  of  average  data  was  identified  bv  a 
sequence  number  and  a  calendar  date.  These  data  were  printed  and  also 
stored  in  a  sequential  data  set  on  a  magnetic  disk  pack.  The  disk  data  set 
could  then  be  accessed  by  SUPERMOCK  to  obtain  main-stem  river  stage  everv  10 
days  for  the  duration  of  a  simulation  period. 

Preconstruction  and  postconstruction  main-stem  data  sets  were  created 
by  RIVCHANGE.  Differences  in  t he  preparation  of  preconstruction-  and  post¬ 
construction-area  data  sets  involved  accessing  different  master  data  sets 
and  specifying  a  dilterent  set  of  associated  grid  nodes  and  river  miles. 

Time-variant  10-dav-average  stages  on  significant  tributaries  to  the 
Red  River  were  also  required  by  SUPERMOCK.  A  digital  program  called  TR 1 BCHANGH 
was  developed  to  provide  these  data  in  a  suitable  form.  Input  requirements 
for  TRIBCHANC.E  include  the  following:  (1)  the  total  number  of  tributary- 
stream  nodes  to  which  stages  would  be  assigned,  (2)  a  beginning  sequence 
number — identical  with  that  for  the  main-stem  data  set-for  computing  sequence 
numbers  for  sets  of  tributary-stream  output,  (3)  manually  computed  10-dav- 
average  stages  at  gaging  stations  on  each  stream,  and  (4)  associated  grid 
nodes  and  stream  miles  for  eacli  stream.  Data  for  any  number  of  streams  can 
be  used  as  input  to  TRIBCHANGE,  and  the  entire  tr ibutarv-stream  data  set  may 
be  created  in  one  run  of  the  program. 

TRIBCHANGE  was  designed  to  interpolate  onlv  spatially  because  tiie  10-dav 
averages  entered  to  it  were  computed  manually  for  the  needed  time  increments. 

At  nodes  where  tributary  streams  enter  the  Red  River,  the  10-dav-average  data 
from  tiie  main-stem-data  set  were  entered  to  TRIBCHANGE  as  data  for  the  base 
gage  on  that  stream. 

Output  from  TRIBCHANGE  consisted  of  10-dav-average  stages  everv  10  davs  at 
all  specified  grid  nodes  for  tributary  streams  in  a  particular  lock-and-dam 
area.  Each  set  was  identified  by  a  sequence  number  identical  with  the  sequence 
number  of  a  corresponding  average  set  in  the  main-stem-data  set. 

Data  from  TRIBCHANGE  were  printed  and  also  stored  in  a  sequential  data 
set  on  a  magnetic  disk  pack.  This  disk  data  set  was  then  accessed  bv  Sl'PKRMOCK 
to  obtain  10-day  average  tributary-stream  stages  everv  10  davs  during  the 
duration  of  a  run. 

Both  preconstruction  and  post cons t rue t ion  iributarv-stream-data  sets 
were  created  by  TRIBCHANGE.  Changes  in  data  used  as  input  to  the  program 
for  postconstruction  included  changes  in  base-gage  data  at  the  mouth  of 
streams  emptying  directlv  into  the  Reel  River,  thereby  reflecting  increased 
postconstruction  stages  on  the  Red  River. 
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in  an  a  1 1’iiame  r  i  e  nap  ol  the  grid.  The  a  I  pit. utter  i  e  nap  tor  I  r '.Ml  was  j  dent  ieal 
with  tile  map  used  in  SITEKNOCK  tor  uletit  i  I  v  ini’,  ev.ipot  r.msp  i  rat  i  on  and.  h'.airanlie 
1  'oiidne  t  i  v  i  t  v  ot  the  eonl  ining  lavi-r.  (See  tie..  It-.  ) 

Values  ol'  AET/AH  were  de  t  e  nn  i  tied  with  the  aid  ot  a  d  i  a  i  t  a  I -eonpu  t  e  r 
program  e. tiled  DKl.KTHKI.H.  Input  data  requirements,  program  list  in,:,  and  an 
example  ol  program  output  are  shown  in  at  t  aehnient  K.  rite  eomputation  seheme 
in  this  program  is  based  on  a  method  given  hv  Ripple,  Kuhin,  and  van  livlekama 
(1S7J).  I'rimarv  input  to  the  program  includes  tlu-  following:  (1)  the  uppi  r 
it  nil  lower  hvdraulie  eondtie  t  i  v  i  t  i  es  OUT  and  HCL)  ot  tin-  con  f  i  n  i  tie,  laver  in 
each  subarea,  as  determined  from  calibration  of  Sl'PEKMOt'K;  (J)  the  thickness 
ot  material  from  tlu-  base  ot  the  root  zone  to  tin-  top  ot'  tin-  aquifer  in  each 
subarea;  and  (1)  values  of  evapotranspi rat  ion  divided  hv  saturated  hydraulic 
coiuluc  t  i  v  i  t  v  (EWETO)  for  depths  to  the  zone  of  saturation  of  from  1  to  it)  ft 
(0.  (  to  4.1  ni)  lor  four  ranges  in  hvdraulie  conductivity,  as  computed  hv 
A  TMOIITK .  ising  these  data,  DE1.ETDEI.H  computes,  lor  eacii  subarea,  values  of 
AET/AH  tor  depths  to  water  of  from  l  to  30  it  (0.  1  to  d.l  m)  . 

Values  ol  AKT/AII  are  computed  l>v  hEI.ETDELH  hv  the  following  procedure. 

Input  valtu-s  of  (1WET0  (ratio  of  limiting  rate  of  evapotranspi rat  ion  to  hvdraulie 
conductivity)  were  multiplied  by  the  input  value  of  1KT  to  convert  tlu*  dimen¬ 
sionless  (1WKT0  value  into  a  limiting  rate  ot  i-vapot  ransp  i  rat  i  on  (in  feet  per 
day)  from  tin*  water  table.  This  flow  was  then  routed  down  to  the  base  ol  tin- 
confining  laver,  using  input  data  on  MCI.  and  thickness  of  the  confining  laver, 
to  obtain  tlu*  artesian  head  (expressed  as  depth  to  water)  neeessarv  to  sustain 
this  I  Low.  The  flow  is  largest  for  the  shallowest  computed  depth  to  water  (1 
ft,  or  0.3  m)  and  decreases  with  increasing  depth  to  water.  The  steadv-state 
model  uses  changes  in  head  and  flow  as  boundary  conditions.  Change  in  flow 
per  unit  of  head  change  (AET/AH)  was  computed  hv  DEEETDEI.H  for  an  input  value 
ot  deptli  to  water  hv  dividing  differences  in  flow  by  d i f f erences  in  depth  to 
water.  An  example  of  results  of  the  computations  is  shown  in  attachment  T. 

Tlu*  relation  between  evapot  ransp  i  rat  ion  and  deptli  to  water  is  .1  curvil  i- 
noar  function.  The  function  is  computed  hv  the  program  DELETPEEH.  Output 
from  this  program  are  tables  of  AET/AH  values  for  depths  ot  I roni  I  to  10  it 
(0.  1  to  l) . L  m)  to  water.  The  model  calculates  the  change  in  evapot ransp i ra t ion 
witii  change  in  water  level  as  a  linear  t unction.  Therefore,  an  iterative 
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1 1  vas  assured  that  the  stages  of  tributary  streams  will  regain  anchor.- 
hrouah  t:.e  cost  con.- tract  ion  period,  except  where  affected  by  backwater  from 
avication  pools.  At  t:.e  r.ioutii  of  each  tr ibutarv  stream,  tine  chance  in 
ate r-sur face  elevation  was  set  ecual  to  the  chance  in  the  rain  sti-i  at  t:;at 
oint .  I'n  i .-  chance  in  water-surf  ace  elevation  was  extrapoiat-d  upstream  to 
oint-of -zero  chance  win-re  the  protected  backwater  profile  intersected  the 
atural-st  rear,  prof  tie. 

cAi.ihKAnoy  AN.)  V.-.RI.-  1CA1  ION  i!i  i  Hi.  NON  S  Tr.A.)!  -  STA  .  1.  Mi'.N'.I. 

Calibration  and  verification  (V.atalas  and  'haddock,  l^T*)  of  the  nonstea 
cate  model  wa-  'cased  or.  simulating  the  observed  water-table  levels  and 
otent  ior.etr  io  heads  at  observation-well  locations.  Calibration  was  ef  fecte 
iti. in  esta:- 1  i  sited  error  criteria  by  adjusting  rodel  parameters  wititin  estab 
i-:.eu  plaus ibi 1 ity  ranees.  Computations  were  made  usinc  river-stace  and 
lir:at  io  data  for  a  -—.n-ar  period  of  record  to  provide  sufficient  tire  for 
nclu- ion  o:  antecedent  conditions.  Visual  inspection  of  computed  --wear 
i  ro  ;  rap :.  s  indicated  that  most  of  tile  antecedent  conditions  were  satis;  ied 
urine  tin/  first  calibration  year  and  that  all  had  been  satisfied  by  the  end 
f  t::/  second  war.  line  tit  i  rd  and  :  our  tit  years  o;  the  observed  water-level 
a t a  were  sniit  into  two  periods.  Tine  fourth  year  was  chosen  as  the  calihra 
; .  >n  juried  for  obtaining  a  rate!:  of  simulated  and  measured  ir/drocraniis  hv 
;  .-tin.:  model  parameters.  The  rodel  results  tor  the  third  'ear  were  used 
.  r  /■  r  i  f  i  cat  ion  evaluation  bv  compar  i  niz  dill  erences  between  t:te  measured  an 

;  .1  it*  i  h  "d  ro  c  ran  :  s  . 


In  the  calibration  process,  all  computed  water-table  and  potent iormtric- 
surtaee  elevations  ot  observation-well  locations  in  the  grid  were  passed  bv 
SUPERMOCK  on  a  magnetic-disk  data  set  to  DATE.  DATE  performed  several  func¬ 
tions,  including  (in  sequence)-'  (!)  converting  mean  sea  level  elevations  to 
depths  below  land  surface;  (2)  assigning  calendar  dates  to  .ill  water  levels; 
(5)  choosing  the  spring  high  and  the  tall  low  water  level  for  tlu-  water  table 
and  potent iometric  surface  for  1  or  more  Years,  as  specified;  (4)  comparing 
these  high  and  low  computed  values  with  observed  data  entered  to  it  on  cards; 
(5)  printing  a  calibration  table  for  analysis;  and  ((>)  passing  all  computed 
water-table  and  potent iometr ic-surface  levels  lor  the  observation  nodes  in 
card  images  to  HYDROG  (Reed  and  others,  197b)  on  a  magnetic-disk  data  set.  An 
example  ot  the  calibration  table  produced  by  DATE  is  shown  in  figure  20. 

Using  these  computed  data  as  input,  HYDROC  plotted  hvd rographs  for  both  the 
potent iomet r ic  surface  and  the  water  table.  Those  computed  hvdrographs  were 
compared  visually  with  observed  hvdrographs  to  check  for  differences  between 
the  two  in  fluctuations  and  depths  to  water. 

The  parameters  that  were  modified,  within  predetermined  plausibility 
ranges,  during  calibration  of  the  nonsteady-state  models  included  the  upper 
and  lower  hydraulic  conductivities  of  the  confining  layer  (HCU  and  HUE)  in 
each  subarea  of  the  grid,  aquifer-storage  coefficient  in  each  subarea  (S), 
specific  yield  in  each  subarea  (WTSTO),  and  streambed  thicknesses  (AM).  In 
order  to  match  observed  data,  a  general  calibration  table  from  DATE  was 
inspected  to  determine  in  which  subareas  simulated  water-table  and  (or) 
potent iometr ic-surlace  levels  were  within  responding  predetermined  error 
criteria.  In  addition,  computed  potent iometrie-surfnee  and  water-table 
hvdrographs  from  HYDRtXi  were  compared  to  hvdrographs  developed  from  observed 
measurements.  After  a  thorough  analysis  of  a  run,  indicated  changes  were  made 
to  appropriate  parameters,  and  a  new  computer  run  was  made.  Normal lv,  from  2(1 
to  25  runs  were  required  to  calibrate  each  nonsteadv-state  model. 

The  magnitude  and  direction  of  changes  that  were  caused  by  modification 
of  parameters  during  calibration  of  the  nonsteadv-state  model  are  discussed  in 
the  following  paragraphs. 

Tile  hydraulic  conductivity  of  the  upper  segment  of  the  confining  layer 
(HCU)  limits  recharge  to  the  water  table  and  thus  the  aquifer.  Eor  positive 
accretion,  an  increase  in  the  modeled  value  of  HCU  in  a  particular  subarea 
causes  an  increase  in  elevation  of  the  water  table  and  very  likclv  an  increase 
in  the  elevation  of  the  po ten t iomet r i c  surface  unless  the  conductivity  of  the 
lower  segment  of  t lie  cent  ining  laver  (IICL)  is  very  low.  A  decrease  in  the 
modeled  value  of  HCU  has  the  opposite  effect  and  would  likely  cause  a  decrease 
in  the  elevation  of  both  the  water  table  and  the  potent iomet ric  surface. 

As  previouslv  mentioned,  HUE  is  the  designation  ot  tile  hvdratilic  conduc¬ 
tivity  of  the  confining  layer  from  the  water  table  to  the  top  of  the  aquifer. 
An  increase  in  the  modeled  value  of  HUE  general !v  results  in  an  increase  in 
the  elevation  of  tiie  po ten t i ome l r i c  surface  and  a  decrease  in  the  water  table. 
Decreasing  the  modeled  value  ol  HUE  lias  the  opposite  el  feet;  that  is,  the 
water  table  will  rise  and  the  po Lent  i ome t r  ic  surface  will  tall. 
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It  negative  accretion  is  occurring — that  is,  the  water  table  lias  dropped 
below  the  potent iometric  surface  and  water  is  moving  up  from  the  aquifer — 
increases  in  the  modeled  value  of  HCL  will  raise  the  water  table  and  lower  the 
potent iometr ic  surface.  Decreases  in  the  modeled  value  of  HCL  will  lower  the 
water  table  and  raise  the  potentiometri c  surface.  Changes  in  the  modeled 
value  of  HCL'  will  either  raise  or  lower  both  surfaces. 

Adjustments  to  the  storage-coef f icient  (S)  and  speci f ic-vield  (WTSTO) 
values  can  be  used  to  control  the  fluctuations  of  the  potent iometr ic  surface 
and  water  table,  respectively.  Increases  in  modeled  values  of  either  of  the 
two  parameters  will  cause  smaller  fluctuations,  and  decreases  in  these  values 
will  cause  larger  fluctuations.  Modifications  to  these  narameters  are  very 
useful  during  calibration  for  adjusting  computed  water  levels  for  the  spring 
and  fall  that  differ  from  observed  values  by  about  the  same  magnitude  but  in 
opposite  directions. 

The  thickness  of  streambed  materials  can  be  adjusted  to  produce  a  change 
in  the  comnuted  water-table  and  potent  iometr ic  levels  near  the  stream.  In  the 
Red  River  Valley,  the  movement  of  water  generally  is  from  the  alluvial  aquifer 
to  the  river  or  stream.  Therefore,  if  the  modeled  thickness  of  streambed 
material  is  too  small,  computed  water  levels  in  observation  wells  near  the 
stream  are  lower  than  the  observed  levels  unless  the  recharge  rate  is  in¬ 
creased  substantially.  If  increasing  the  conductivities  of  the  upper  confining 
layer  in  the  affected  subareas  does  nothing  more  than  increase  computed  accretion 
values  to  unreasonable  levels  without  an  appreciable  rise  in  water  levels,  the 
streambed  thickness  is  too  small.  The  upper  plausibility  limit  on  accretion 
was  1  ft/yr  (0.3  m  per  year).  An  example  of  a  table  of  annual  accretion 
summations  at  observation  wells  is  shown  in  t  igure  21. 


MODEL  Ol'TPl’T 

The  output  from  the  nonsteadv-state  model,  in  addition  to  that  used  for 
calibration,  was  designed  to  display  the  results  of  the  analysis  in  a  form 
suitable  for  the  determination  of  the  effects  of  the  water  table  on  agricul¬ 
ture.  The  critical  parameter  influencing  agricultural  production  is  the  depth 
to  the  water  table  below  land  surface.  The  depth  to  tin-  water  table  has  a 
significant  effect  when  it  Is  within  the  root  zone,  or  within  approximately  3 
ft  (1.5  m)  of  the  land  surface.  Times  of  occurrence  of  shallow  depths  to  tin- 
water  table  are  also  significant.  The  most  critical  periods  occur  during,  tin- 
plowing,  planting,  growing,  and  harvesting  seasons  of  the  vear.  for  this 
reason,  output  from  the  model  was  designed  to  show  the  average  depth  to  the 
water  table  for  either  one  or  two  10-dav  time  frames  during  these  critical 
periods.  A  series  of  30-day  time  frames  was  used  to  represent  water-table 
conditions  during  the  dormant  season.  In  this  manner,  the  vear  was  divided 
into  21  time  frames  associated  witii  specific  calendar  dates.  The  dates  were 
selected  by  the  Soil  Conservation  Service.  Tin-  actual  output  consisted  nt 
data,  punched  on  computer  cards,  showing  the  computed  depth  to  the  water  table 
below  land  surface,  to  the  nearest  foot,  at  each  node  in  the  model.  figure  22 
is  an  example  of  part  of  the  data,  in  printout  form,  showing  the  node  local  ion 
(row  and  column)  and  depth  to  the  water  table  for  a  particular  tine  trame. 
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Figure  22. --Example  of  computed  output  from  nonsteady-state  model. 


Depths  to  the  water  table  of  as  much  as  9  ft  (2.7  m)  below  land  surface  are 
shown  in  the  printout.  Because  the  table  format  prints  depth  as  only  a  single¬ 
digit  integer,  a  9  indicates  a  depth  of  9  ft  (2.7  m)  or  more. 

The  combination  of  depth  to  water  table,  soil  type,  and  cropping  pattern 
was  used  by  the  Soil  Conservation  Service  to  determine  the  beneficial  or 
adverse  effects  of  project-induced  changes  in  water  levels.  Crop  yields 
obtained  during  the  calibration  period  were  used  as  the  standard  for  determining 
the  net  project  effects. 

Output  from  the  steady-state  analysis  consisted  primarily  of  head-change 
data,  shown  as  tabulations  or  as  maps.  Maps  of  head  change  with  time  are 
available  from  the  model,  but  only  the  final  or  steady-state  output  was  consid¬ 
ered  significant  because  it  represented  the  dynamic  equilibrium  conditions 
resulting  from  the  change  in  river  stage.  This  head-change  map  was  used  to 
compute  the  average  postconstruction  potent  iometric  surface.  The  elements  of 
this  computation  are  shown  in  figures  23-26.  Figure  23  shows  the  preconstruction 
po t en t  iome t r  i c  surface  in  a  lock-and-dam  area.  The  computed  head  change  is 
shown  as  a  grid  plot  to  the  same  scale  as  the  model  grid  (fig.  24)  and  as  a 
complete  contour  map  (fig.  25).  The  head  change  was  added  algebraically  to 
the  preconstruction  potent iometric  surface  to  produce  the  resultant  potentio- 
metric  configuration  shown  in  figure  26.  This  method  is  based  on  the  principle 
of  superposition  that  assumes  that  the  flow  field  in  the  aquifer  can  be 
considered  a  linear  system  and  that  the  head  change  component  can  be  analyzed 
independently.  The  principle  of  superposition  allows  the  postconstruction 
condition  to  be  determined  as  the  sum  of  the  preconstruction  head  and  the 
head-change  component. 


CONTINUING  STUDIES 

The  modeling  procedures  developed  for  this  study,  particularly  those  for 
the  modeling  of  nonsteady  flow,  were  designed  to  provide  data  lor  an  assessment 
of  the  effects  of  proj ec t- induced  water-level  changes  on  agriculture'.  However, 
these  procedures  can  be  applied  to  a  variety  of  situations  in  connection  with 
the  Red  River  Waterways  Study.  The  calibrated  models  can,  with  the  appropriate 
boundary  changes,  be  used  to  analyze  the  effects  of  any  arrangement  of  locks 
and  dams  or  pool  stages.  Although  the  results  of  the  study  wore  primarily 
concerned  with  agriculture,  the  nonsteadv-state  model  can  be  modified  to 
determine  the  effects  of  raised  water  levels  in  urban  areas.  The  higher  water 
levels  may  cause  flooding  of  basements,  septic  tanks,  or  sewer  systems,  or 
may,  because  of  increased  moisture'  content  of  surficial  clays,  cause  differen¬ 
tial  movement  of  footings  of  buildings,  swimming  pools,  or  bridges.  The 
models  can  also  be  used  to  aid  in  the  design  of  well  fields  and  surface- 
drainage  systems  that  mav  be  needed  in  places  where  shallow  water-table 
levels  are  anticipated. 

To  achieve  the  greatest  benefit  from  the  studv,  the  water — love  1 -obsorva t i on 
network  developed  for  the  studv  should  be  maintained  and  water-level  measurements 
continued  through  the  construction  phase  to  verify  the  predictions  made  during 
the  study.  The  data  would  provide  a  definition  of  the  actual  ground-water 
conditions  resulting  from  the  stage  changes  and  would  provide  a  means  of 
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0 .8 

0. 1 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

3.1 

2.9 

2.7 

2  •  J 

2.1 

1.9 

1.8 

1.3 

0.6 

0.2 

0.  1 

0  .  C 

0 . 0 

0 . 0 
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0 . 0 
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1.0 
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0.0 

0.0 
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0.0 

9.5 
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0.7 

0.5 

0.4 
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0.0 

0.0 

0.0 

0 . 0 

0 . 0 

0.0 
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5 .5 

3.3 
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0.0 

0.0 

0 . 0 

0.0 

0.0 

0 .  U 
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4 . 1 
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0.3 

0.2 
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0.0 

0.0 

0.0 

0  .  u 

0.0 
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0.0 

0 . 0 

vJ  •  u 
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2.4 
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0.1 
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0 . 0 

0 . 0 

0 . 0 
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NOTE. —  Eight  sheets  of  printout  required  for  complete  covera'i"  of  .1  sinrlc 
modeled  area. 
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ATTACHMENTS 

fhe  following  attachments  give  the  program  listings  and  show  the  input 
data  requirements  for  the  peripheral  programs  used  in  conjunction  with  the 
Sl'PEKMOCK.  and  GWFI.OW  models.  The  relation  of  the  peripheral  programs  to  the 
models  in  shown  in  figure  27.  Examples  of  printed  output  (figs.  28,  29)  from 
the  AVERAGE  and  DKLKTDELII  programs  are  given  along  with  the  respective  docu¬ 
mentation.  Output  from  the  AVERAGE  program  is  not  used  directly  in  the  GWFI.OW 
model,  but  it  provides  control  for  tile  contour  map  of  the  average  preconst ruc- 
t ion  potent iometrie  surface  used  in  conjunction  with  the  output  from  the 
GWKLOW  model.  Primarv  output  from  the  ATMOFLl'X  and  pnTRKT  programs  is  punched 
on  cards  and  that  I  rum  the  R1YCHANGF.  and  TRIBCHAN'CE  programs  is  stored  on  disk 
data  sets.  Examples  of  output  from  these  four  programs  are  not  shown. 
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Card  input  containing  location  I  I  faster  main- 


ATTACHMENT  A 
AVERAGE  Program 
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Table  3. 
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A  V  F  ►-  A  S  F 

p  Y 
Tv 


THIS  PPOOPA*  CO^PUTF  S  T  “y  mII'Jy  0  r-t  I r  [TlYi  -.t  tAhTi-i  i.  a  of  ,A 

L  K  V  F  L  FOP  A  SPFCIFIFD  PFPJTlI.  OF  PFO>  P'‘  Ri,-  ,■ .  ,  v  y  '  I  y  >  .v  v  .nr,, 
Thf  ppoopam  is  oRM'vNFp.  to  Yfin  |  flY',  *,  p  rA.."s  is  r-t  .,  <. 
GPOiiND-pATFp  F0PMATS.  Hfu»F\/FP.  IT  ro.||'.  pc  ..ff  .|H<aity  t  < .  :.((••• 

WATFP-IFVFL  CftPOS  AS  Tuf  V  A  P  F  MlINlfpC  l  F  f,„  v  y  S  T  F 

F  in>  (J  o  -J  iHHUV  5  Y  Y  rf  J  Y  O  'J  V  ft  -.t  Y  Y  Y  1Y  Y  Y  f>  •;>  t-  ■.,  :v  .1  >,  ,  -.>  i  .  ....  ...  ^  „v  .,»  ....  .;,.  j  . ,  .,  .  <■  V  .:  v> 

0  I  MF NS  I  ON  I  OA  V  (  1  ?  )  ,S  IOM  (  T  )  ,  I  HFD  (  pno  |  ,  pr  (  a  )  .  !  '■  (  A  )  .  :  y  (  4  !  .  v|  (  u  )  . 

!  MN(4).IMn<4).Iln<4),MY(4).  AVI.  (  4  ) 

DATA  JM, jn.  IV. W|  M/0.0. 0.0.1/ 

DATA  M(,N/'  •  / 

DATA  T0AY/T1  »  ?P  «  3  1  <  in.  01.  3  ■  1.0].  II.  TO. 0|.  I.,.  <  ]  / 

I  P0  =  S 
I  P  T  =  P 
T  7  T  =  0 

I^ON  --  WF'GINNING  «0NTh  FOh  ofsJpF'i  prj;ni|  OF  pfti.ui.. 

I  Y  P  --  HFGINNING  YFAP  FOP  Of  S  I  PF  O  OF  P  I  no  OF  PFOC.PO. 

JMpN  --  f  MO  I  NO  u0NTh  FOP  OF  S  I  pF  o  pFp[oo  of  .trns), 
jyp  __  FajOInG  YFAP  FOP  O  F  S  I P  F  o  Ff  p[ril  OF  PFOOPI.. 

PF An  (IPD.4S)  T«0N. I YP, JMON, JYP 

FlF V  --  L ANO-SUPF Af F  FL  FVAT  TON  AT  aF^I  . 

I M  n  --  CAPO  MJMPFP. 

I  PF  A  n  (  I  P  0 , 4  f,  ,  F  MO  :  4  4  )  Fl.FV.INO 
Mijrfl 

I  FMrfl 
js«  =  n 

DFPTm  --  OFPTH  OF  WFLI. 

IMO  --  CAPO  N'JMPFP. 
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pfao  (TPn.47)  ofpth. inn 
IF  (OFPTH. FO.n.n)  OF PTPrqgp.gp 
JM  =  Ti 


n  o  r»  o  n  n  r>  o 


Table  3. — 41'^''  ;.r.-jri -  ;Ur:>.  .-—Continued 

ww  =  n .  n 

pa  x  =9999  . 

NPr  1 

NTM=n 

NQ  =  1  4 

P!N=-999.9 
A VFl =  999 . 99 
PAXFL=999.99 
PINFL=999.99 

r 

C  I HFO  (  I  )  --  HF  AO  I  NO  INFORMATION  for  wF l  L  • 

r 

C  I N  0  --  CAPO  NIJMHFO. 

c 

?  PFAO  (  I  P  0  *  4  M  )  {  IHFOI  I)  ,I=NP,N'J)  .  INO 

IF  (INO.LT.l)  fiO  TO  3 
NP=NP. 1 4 
NQ=NQ ♦ 1 4 
GO  TO  ? 

3  MQ=MQ*1 

M=] 

c 

“0(1)  --  MONTH  WATFP-LFVFL  “FASUPMFNjT  T4k) N, 

10(11  --  OAY  wATFP-lFvt'L  MFASUPFMFNT  T«OM. 

I  Y  (  I  )  —  YFAf,  WATFP-LFVFL  MEASUREMENT  T  a  *  F  N  . 

S I N ( I )  —  SIGN  OF  WATFP-LFVFL  VALUE. 

r  WL ( J  >  --  WATFP  LF VFL 

n 

4  RF  an  (IPO. 49)  (“0(  I  )  ,  m  (  I  )  .  I  Y  (  I  )  ,SIN  (  I  )  .WL  (  I  )  ,  I  =  1  ,4  )  ,  IMIJM 

IF  (MO. ME  .  I  )  GO  TO  9 

IF  (IYR.FO.O)  GO  TO  9 
00  7  I T= 1  .4 
IF  (IY(IT)-IYW)  7,4,4 
c  IF  (MO ( I T ) .LT . I “ON)  GO  TO  7 
4  IF  (ISM.FQ.l)  GO  TO  7 
I  SM=  1 
M=  I  T 

7  CONTIMJF 

IF  (ISM.FQ.l)  GO  TO  P 
IF  (JM)M.FQ.O)  GO  TO  1 
GO  TO  4 
P  L  MrMO ( M ) 

LOr 10 (M) 

L  Y= I Y (“) 

LPMrO 

9  00  10  I - M . 4 


Alt  4  ] 

All  4  p 

Alt  49 

All  44 

A  V  f  4  4 

A  V  F  ^4 

Alt  4  7 

Alt  4  p 

A  IF  4g 

AVt  41) 

Alt  4  ] 

Alt  f; 

AVI.  4  1 

Alt  44 

Alt  44 

A  l  F'  4  4 

A  l  F  4/ 

Alt  44 

A  IF  40 

Alt  7f| 

AVt  7  1 

AVt  7? 

A  V  K  7  1 

Alt  74 

Alt  74 

Alt  74 

AVt  77 

AVt  74 

AVt  79 

AVt  BO 

Alt  P) 

Alt  RP 

AVt  P  3 

A  Vt  R4 

AVt  H4 

AVt  R  4 

AVt  P  7 

AVt  PP 

AVt  P  9 

AVt.  911 

AVt  9  1 

Alt  9  p 

AVt  99 

AVt  94 

AVt  94 

AVt  94 

AVt  9  7 

AVt  9R 

A  Vt  99 

Alt  100 
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Table  3.—  ’.'-  ?»•<•>  •>*»»•  • — Continued 


IF  (SIN( I ) .FQ.SIGN(  1 > )  WL < I > =WL ( I ) * (-1 .0) 

10  CONTINUE 

IF  (JYR.FQ.O)  GO  TO  14 
I F  ( JSM.FQ. 1  I  GO  TO  37 
po  n  JT=M,4 

IF  (IY(JT)-JYR)  13,11,1? 

11  IF  <MO< JT) .LE.JMON)  GO  TO  13 
1?  JSM=1 

IF  (JT.FQ.M1  GO  TO  37 
MO ( JT ) =MO ( JT-1 ) 
in ( JT 1 =10 ( JT-1  ) 

I Y ( JT ) = I Y (JT-1 ) 

WL ( JT) =  WL ( JT-1 ) 

13  CONTINUE 

14  IF  (M.FQ.4)  GO  TO  36 
PO  33  I =M , 3 

IKO  =  0 
JT  =  o 
I  JJ=  0 
JJJ=0 

IF  (JM.FO.O)  GO  TO  IS 
GO  TO  ?0 

1R  IF  ( MO (1*1), NF . 0 )  GO  TO  16 
MO ( I ♦ 1 ) =MO ( I ) 

IO(I*l)=IO(I) 

IY ( 1*1 )  =  IY  (  I  ) 

WL  (  I  »  1> =  WL (  I  ) 

GO  TO  30 

16  IF  (MO ( I ) .NF .MO ( I ♦ 1 ) .OR. I Y ( T ) .NF  .  I  Y  ( I  ♦  1  )  )  GO  TO  17 
P4=(l,n*(I0(I*l )-IO(I)  )  )/?.0 

WW=WW*  (WL  (  I  )  *0  A )  ♦  (WL  (  1*1  )  *<BA-1 .0)  ) 

NTM  =  NTM*  (10(1*1  ) - 1 H  (  1 1  ) 

GO  TO  30 

17  IF  ( I Y ( I ) .NF.  I  Y (  I  ♦ 1  )  )  GO  TO  ?0 

IF  (MOOClY(I) ,4) . FQ . 0 )  IOAY (?) =  ?0 
K=MO (1*1) -MO ( I ) 

IF  (K.FQ.l)  GO  TO  I R 
I  a  =mo ( T ) *1 
I R=MO ( I ♦ 1 ) - 1 
00  10  J=I A, IB 
10  JJJ=JJJ* IOAY ( J) 

1 Q  NO Y  = I  0 A  Y (MO ( I )  1 -10  (  I  ) 

BA= ( l . o*NOY* JJJ* 10 ( I* 1 ) ) /?. 0 
WWsWW* (WL ( I  1 *BA ) ♦ (WL ( I *1 ) • (BA-1 .0 ) > 

NTM=NTM* JJJ* 10 ( r* 1 ) *N0Y 

GO  to  30 

?0  IF  (  ( IY ( 1*1 ) -I Y (  I)  )  .GT. I )  GO  TO  ?  6 
?1  KRsMO ( I ) * l 

IF  (MOOdY(I)  ,4).FO.O)  IOAY  (?)  =  ?0 
IF  (KP.GT.l?)  GO  TO  ?3 


A VF  101 
AVE  10? 
A  VF  103 
AVE  104 
AVE  10S 
AVE  106 
AVE  107 
AVE  100 
AVE  100 
AVE  110 
AVE  111 
AVE  11? 
AVE  113 
AVE  114 
AVE  ns 
AVE  116 
AVE  1  1  7 
AVE  11« 
AVE  119 
AVE  1?0 
A VF  1?1 
A  VF  1?? 
AVE  1 ?  3 
A VF  1?4 
AVE  1?S 
AVE  1?6 
AVE  1?7 
AVE  1?6 
AVE  I?9 
AVE  130 
AVE  131 
AVE  13? 
A  VF  133 
AVE  134 
AVE  13S 
AVE  136 
A  VF  137 
A VF  130 
AVE  130 
AVE  140 
AVE  141 
AVE  14? 
AVE  143 
AVE  144 
A  VF  14S 
AVE  146 
AVE  147 
A  VF  140 
AVE  140 
AVE  ISO 


Table  3. — AVERAGE  program  listing  —  Continued 


00  2?  J=KP,12 

22  JJJ=JJJ*IDAY ( J) 

IOAY (?) =28 

23  IF  (MO(I*l).EQ.l)  GO  TO  25 
I PP  =  MO ( I ♦ 1 ) - 1 

IF  ( MOD  (IY  (!♦!)  *4)  •  E  Q  •  0  )  I0AY<2)=?9 

00  24  J=1«IPP 

24  I  J J= I JJ* I  DA Y  <  J ) 

25  NO Y= ID AY (MO(I) ) - ID  1 1 ) 

BA= ( 1 . 0*JT*NDY*JJJ*I JJ*ID ( 1*111/2.0 
W  W  =  W  W  * (WL < I )  *8 A ) ♦ ( WL ( I ♦  1 ) * (BA-1 .0) ) 
NTM  =  NTM*N0Y* JJJ* I JJ* ID  <  I  ♦  I  )  ♦  JT 
GO  to  30 

25  KP=(lY(I*l)-lY(I) )-l 
KK I =  365 
DO  27  J=1.KP 

IF  (MOO ( < IY ( I ) ♦ J) .4) .EQ.O)  KK I =366 

27  JT  =  JT *KK I 
GO  TO  ?1 

28  IK0=1 

00  29  K=  1  * 4 
IMO(K) =MO(K) 

IID(K)=ID(K) 

I  I Y  <K ) =1 Y (K) 

AWL (K) =WL (K) 

IF  (K.FQ.l)  GO  TO  ?9 
MO(K) =M0(K-1 ) 

ID (K) =ID (K-l ) 

IY (K) =IY (K-l ) 

WL IX) =WL (K-l ) 

29  CONTINUF 
M0( 1 ) =JM 
IP(])=JD 
I  Y  (  1  ) =  JY 
WL ( 1 ) =WLM 
GO  TO  15 

30  IOAY (?) =28 

IF  ( IKO.NE. 1 )  GO  TO  32 
DO  31  K  =  1 » 4 
MO (K) =  T  MO ( K ) 
in(K)=IID(K> 

IY(K)=HY(K) 

31  WL  (  K  )  =  A  WL  (  K  ) 

I  KO  =  0 

JM  =  0 

GO  TO  15 

32  LLM=MO (1*1 ) 

LLD= 10(1*1) 

LLY=IY  (1*1 ) 

WPD  =  WL  (1*1 ) 


AVE  151 
AVE  152 
AVE  153 
AVE  154 
AVE  155 
AVE  156 
AVE  157 
AVE  158 
AVE  159 
AVE  160 
AVE  161 
AVE  162 
AVE  163 
AVE  164 
AVE  165 
AVE  166 
AVE  167 
AVE  168 
AVE  169 
AVE  170 
AVE  171 
AVE  172 
AVE  173 
AVE  174 
AVE  175 
AVE  176 
AVE  177 
AVE  178 
AVE  179 
AVE  180 
AVE  1 R 1 
AVE  182 
AVE  183 
AVE  184 
AVE  185 
AVE  186 
AVE  187 
AVE  188 
AVE  189 
AVE  190 
AVE  191 
AVE  192 
AVE  193 
AVE  194 
AVE  195 
AVE  196 
AVE  197 
AVE  198 
AVE  199 
AVE  200 


Table  3. — AVERAGE  program  listing — Continued 


C 


33 

CONTINUE 

AVE 

201 

DO  35  N  =  M»4 

AVE 

202 

P I N  =  AMA  X 1 (PIN.WL  <N) ) 

AVE 

203 

PAX=AMJN1 (PAX.Wl <N> ) 

AVE 

204 

IE  (PIN.NF.WL (N) )  GO  TO  34 

AVE 

205 

MA=MO (N) 

AVE 

206 

M0= ID ( N ) 

AVE 

207 

MY=IY (N) 

AVE 

208 

34 

IE  (PAX.NE.WL (N> )  GO  TO  35 

AVE 

209 

NMA=MO (N) 

AVE 

210 

NMD=IO(N) 

AVE 

21  1 

NMY=IY (N) 

AVE 

212 

35 

CONTINUE 

AVE 

213 

36 

JM=MO ( 4 ) 

AVE 

214 

J0=ID(4> 

AVE 

215 

JY=IY (4) 

AVE 

216 

WLM=WL (4) 

AVE 

217 

37 

IE  (JNUM.NE.O)  GO  TO  3 

AVE 

218 

WW=WW*WPD 

AVE 

219 

ntm=ntm*i 

AVE 

220 

AVE=WW/NTM 

AVE 

221 

IF  (ELEV.EQ.0.0)  GO  TO  38 

AVE 

222 

avel=elev-ave 

AVE 

223 

paxel=elev-pax 

AVE 

224 

pinel=elev-pin 

AVE 

225 

38 

IE  (MOD(IZT,3> .EQ.O)  GO  TO  39 

AVE 

226 

WRITE  (IPT.54) 

AVE 

227 

GO  TO  40 

AVE 

228 

39 

WRITE  ( IPT *50 ) 

AVE 

229 

40 

WRITE  ( IPT.51 ) 

AVE 

230 

WRITE  (IPT, 52)  (IHED(I),I=1, NO ) 

AVE 

231 

WRITE  (IPT, 53)  LM,LD,LY,LLM,LLO,LLY 

AVE 

232 

WRITE  (IPT.54) 

AVE 

233 

WRITE  (IPT, 55) 

AVE 

234 

WRITE  (IPT, 56) 

AVE 

235 

WRITE  (IPT, 55) 

AVE 

236 

WRITE  (IPT.54) 

AVE 

237 

WRITE  ( IPT, 58)  DEPTH, AVE.AVEL, PIN, MA.ND, MY, PAX, NMA, NMD, NMY.PINEL 

.PAVE 

238 

1  AXEL 

AVE 

239 

IE  (PAX)  41,42,42 

AVE 

240 

41 

WRITE  (IPT, 57) 

AVE 

241 

GO  TO  43 

AVE 

24? 

4? 

WRITE  (IPT.54) 

AVE 

243 

43 

WRITE  (IPT, 51) 

AVE 

24* 

WRITE  (IPT, 59) 

AVE 

245 

IZT=IZT*1 

AVE 

246 

GO  TO  1 

AVE 

247 

44 

STOP 

AVE 

248 

AVE 

249 

45 

FORMAT  (412) 

AVE 

250 

7U 
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iqure  ?8 --Fx.imple  ot  output  from  AVERAGE.  progr.i 


ATTACHMENT  B 


ATMOFLUX  Program 


order . 


Table  5. 


CT:P»0 C  OPTIONS  !  [N )  ; 

/it  e  (HHHliHH  ■»  «  »  ilHHHHH  &  i*  i  5  4  t  i  it  £■  iHi  v  v  it  it  it  it  -,t  it  it  -.i  it  it  it  -.t  it  it  -t  it  it  it  >  it  ••••  it  ■ 


ft  Tt-O^L  '  ix 


T  h  I s  PPnoHftti  computes  potential  t  i  P  w  a  » n  wnvf  wr,M  (  nos  t  . 

F  VAPOTPANSp  I -A  T  ION  AT  Top  LAND  FoVFirf  )  rnp  PPUTiS  To  I-,f  *,-T^ 
TAPLF  FP()M  I  TO  50  PT.THf  Pponuji- 

COMPUTES  ThP  HAPmONIC~mFAN  hyQPAULIC  CONDUCTIVITY  fop  I  ftYF“P! 
MftTFPlAL  HY:  KSAT(HAt-.  M  p"  AO  HYp.  COOT)  .  )  =  SUM  (  T  H  [  C  OF  S'*  )  / 
Sl|M(THlCI'NFSS»HYO,CONiJ.)  .  IT  SSSIOOS  A'  A  *  J  p  Sl/P  vAl.IJrS  T0  A 

LOO  OFPFN'DINO  ON  CALCULATED  VALUES  of  liPAT  A-jIT  I  '\|  P  i„  T  V'-UJFV  OF 
FXP.FUFTION.  AN!)  F  Y  P_L  I  M  I  T  .  TH  I  S  p  POOP  AM  l.l  OF  s 
N(GAPQNF»'S  FXPOnF.Nt).  AND  S  1  /  ?  (  T  E  NO  I  ON  AT  *"-ICh 

IINSATUPATEO  H  .  C  . /S  A  T'JRA  TF  0  H.C.  =  1/P)  TO  COMPUTE"  V  F  P  T  (CAL  Fl'wi 

AS  A  FUNCTION  OF  QFPTH.  THE  FUNCTION  I  S  F  (J  .  PT(P'.AO)  ~F 
PIPPI.F,  FT  AL.  WSP  P014-A. 

it  it  it  it  it  *  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  t  ,t  it  it  it  it  it  it  it  it  it  it  it  it  it  it  w 

■>  / 

OCL  HC ( 1 R )  INI T (< 1 5) 0  .)  «F TmF (?: R  ,  TO )  .  L  CO  (  ?<-  )  ,  Th  (  PS  j  , 

F(ICT  ION  (  1  n  )  INI  T  (  (  1  n  )  0  .  )  .FXP_L  IM  I  T  (  c  )  .  (  L  .  *  SAT  )  OFr  P!  A4T(H  , 

F  X  P  (  ]  0  )  RTNAPY  FIXE  0(15,0)  »  «Y  F  L  L  _  N  O  Ch  ah  (  4  )  ,  PUNCn  Flit  output: 

ON  FNnFILP (SYSIN)  GO  TO  PIP 
OFT  F ILF (SYSIN)  EDIT 
I NHC •  ( I  *  HC  (  I )  DO  J=1  TO  NMC)) 

NHC  -  NUMPFO  OF  HYO.  CONP.  TO  HF  np  AP  ( A  X  T  hi  jm  =  ]  c, ) 

I  -  nunrFP  PEFFPP I  NO  TO  A  I  ITHOLOOIC  Typf.  may  8F  «pw;t^«FI|  y 
CHOSEN  WITHIN  THE  P ANGF  OF  1-]S. 

HC  (  I  )  -  HYO.  COND.  (FT/OAY)  OF  l  I  T  hoi.oo  I  r  Typf  I. 

/ 

(COL  (1  )  (P (?)  • S K I P  (  1)  ,  (NhC)  ( F (?)  «F(H,o)  )  ) 

(NEXP,  (FXP  (  I  )  .SUCTION  (  I  )  00  1  =  1  TO  NPXP)  ) 

NFXP  -  All  JMPFP  (  Mfl  X  I  M(JM  =  10)  OF  INPUT  VA|_UFS  OF  N. 

FXP(I)  -  V  AI_  UF  OF  N  CORRESPONDING  TO  SUCTION!  I  )  . 

SUC  T I  ON ( I )  -  VALUE  OF  Sl/P  (In  FEET)  CORRf SPONO  1  NO  TO  exd(I). 

/ 

( COL  (1  )  .F (?)  , S  K I P ( 1  )  ,  (NFXP)  ( F ( P )  . F ( H )  )  ) 

((FXPLIUIT(I)  00  1  =  1  TO  NF  x  P- 1 )  ) 

FXP_LImit(I)  -  UPPPW  LIMIT  OF  KSAT  fop  FxP(I).  mimpm  of  VAIUFf 
IS  NF  XD-  1  (MAX  =  0).  FXP_LIMTT  VA|.:jF'S  M||QT  lit  AP~AUOFO  IN  AOCf'jOI\( 
ORPFR,  (SMALIFST  first  and  L  AP(-P  ot  1  APT  )  .  SINCE  T-iPPP  if 
A  COPRFSPON'PF  NCF  RFTWEFN  F  x  P  {  r  )  A  NO  P  X  R  _l  I  M  I  T  (  I  )  ,  FXF(I)  wIlL 
ALSO  RP  COOFO  IN  ASCPNOINO  ORDER. 

/ 

(COL  (  1  )  .  1  0  F  (R.M  )  ; 

PUT  file (SYSPPINT)  FAIT 

(•H.C.  (.  I  m  I  t  «  .  «  S  1  /  P  •  .  *  F  x  P  0  ME  N  T  •  ,  •  (  F  T  /  "  A  Y  )  *  .  *  (  P  T  >  •  ) 


7(i 


Table  5. — ATMOFL'JX  program  lie-ting  —  Continued 


(  x  (  1  ? )  .  a  .  /  { 4  )  .  A  .  fit  (  I  ;  ,  a  ,  x  (  f  ;  ,  a  .  v  i  4  >  .  a  > 

(  ( t  XP  (  I  )  .  •  <  •  .  f  XF_l.  IW1  M  [  1  .  S'n'T  I  ON  (  I  >  00  1-1  To  -;t-  *  p  _  j  )  5 

(  (  n,f  v  p  _  ]  )  (  r  0 1  (  1  )  .  X  (  3  )  ,  F  l  ?  )  .  X  <  7  >  .  A  .  F  (  M  .  4  )  .  v  (  4  ,  .  t  (  4  )  )  ) 

(  c  X  p  ( Mf;  x  p  )  ,snrr]  Oo  ( mf  >.  p  )  > 

(rot.  ( 1)  .  x  (  3)  .p  ( P)  .  x  (  p  )  .c  (  -  j  )  : 
c  1 1 C  T  I  0  N  =  7  o  .  4  ^  S  0  r  T  I  O  J 
no  k  =  i  to  n  f  x  u  ; 

MFXP(X)  ; 
o l ?  =  SUCT  TON ( k )  : 

C?  s  « 

m=w-i  ; 

c'-3.l41oo/(pNjosiN(3.14lou/p\  )  >  ; 

<  =  S  1  ?-"-F  /TO  .  4  q  ; 
if  x < i .  then  x  =  i .  : 
no  t - ]  m  3 o  ; 

(  =30 .4<M  i  : 

A  =  (  c,  ]  ?  *  r  /  L  )  *•  *  N  : 

no  .1=1  to  ieo: 
x  \  =  x  o  »  n  : 
x  n  i  =  x  n  /  x  ; 
x  m  p  =  x  n  ]  /  x  : 

u-  (  *N-  XN  1  -  A  )  /  (  N*x*  1  -N  ]  X  '.<?  )  : 

x  =  x-u ; 

IF  ,|<0.  TMF  Ki  t!=-U! 

IF  ')<  3  .F  -F  T  HF  ► 1  00  TO  of; 

F_  m  n : 

,J=  1  on  ; 

r? :  f ix f (N, i ) =x-i  .  : 

P  N  n : 

*=  1  s 

L  =  1  o  : 
no  j= i , 3 : 

PUT  filF(PUNCh)  FOIT 
(  (F  RF  (N,  I  )  00  I  =f  ,{  )  ,m) 

( OOL  (l).)0F(7.h).<(4),» m= t ,f  <  p)  )  : 

K=K+ 1 o  ; 

L  =  L  ♦  1  o  ; 

F  MO  ; 
f  Mh : 

A  1  :  OF  T  F  ILF  ( PYSI N)  F  0  I  T 

(WELL  NO.  (|.rn  (  I )  ,  th  (  i  )  on  1  =  1  to  )?n 

/* 

v»fi_l  NO  -  wF  L  L  mumqfp. 

L  C  0  (  I  )  -  L  T  T40L0M  r  -  T  TPF  N  i  I  m  p  f  P  Ffi-’  T-u-  Mh  i'NTT  IKl  Th-  l  m, 

TH  (  I  )  -  THIOk-MF  CF  (FT  )  OF  T  -<F  I  T  h  l  m  I  T  M  TM  |jV,, 

/ 

(  rot  (  ]  )  .  X  (  4  )  .  A  (  4  )  .  ]  ?  (  F  (  P  )  ,  F  {  4  )  )  )  ; 

IF  T h  ( 1  ? )  > o  .  T HF  S'  on  : 

r-FT  FII  F  (SYMN)  F  )  I  t 
(  ( i  rn  <  I  I  .  TH  (  I )  oil  r  =  m  t n  pm  ) 

/«• 

IF  TH  (  1  ?)  >11  TH*-  A  [T  I  3  MF  r  F  s  c-  A  H  v  1-1  H  A  v  F  A  F  -  ffi*  o  (  ,*  F  I  (  F  A  P  F 
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I  is  liny  —  Conti  nuecf 


Table  S.—  ATMOFLUX  i-v-jgvor. 


BLANK)  FOR  LCOtTH. 

(COL (  1  )  »  1  3  (F (?) *F  (4) )  )  * 
p’NO: 


*/ 


THCf THK=0.  S 
00  1=1  TO  255 

THICK=TH (1)5 

IF  TH  I  CK<  =  0 •  THEN  00  TO  A2S 

thk^thk+thick; 

HYD_COND=HC (LCD ( I ) ) 5 
IF  HYDCONO<=0.  THEN  005 

PUT  FILE (SYSPRINT)  EDIT 

(•WELL  NUMBER  '*WELL_NO*«*  UNIT  •  ♦  I  *  '  •  CODE  =  *. 
•,  HYDRAULIC  CONDUCTIVITY  =  O') 

( PAGE  t Af A (4)  iAiF (?)  »A»P<?)  .A)  5 
GO  TO  A 1 5 
ENOS 

THC=THC*THICK/HYD_CONOS 
END  5 

A?:KSAT=THK/THCS 

no  i=i  to  nexp-15 

IF  KSAT<EXP_LIMIT ( I )  THEN  DOS 
N  =  F  XP ( I )  S 
GO  TO  CIS 
ENOS 

FNDS 

N=EXP (NFXP) S 

Cl : PUT  FTLE (SYSPRINT)  EDIT 

(•WELL  NIJMRER  •  t  WELL  _N0  ,  •  S  A  TUP  ATED  HYDRAULIC  CONDUCTIVITY 
•  FT/DAY* , 'GARDNER* *S  EXPONENT  =  «,N) 
(PAGF,A,A(4),C0L(1)»A,F(7*4),A,C0L(1)*A.F(?)) 

( 'DEPTH* , »T0» « 'WATEP' (FT) •  ET/SHC  •»* FT ( FT/DA Y )' ) 
(SKIP  (?)  »  A  »  COL (1)  iX(?)  ♦  A  »  X  (  1  B  )  *  COL <1>*A»X(?0>» 

COL  <1)«X(1)»A#XU)»A«X(4).A) 

( ( I .EINF (N, I ) ,KSAT*ETNF (N, I )  DO  1=1  TO  30)) 

(30  (COL  (1  )  »  X  (  2  )  »  F  (?)  » X  (S)  iF  (  B  *  E> )  »  X  (  F  )  »  F  (8*5)  )  )  S 

PUT  FILE(PUNCH)  EDIT 

(WELL  ..NO,KSAT,KSAT,THK,  »  .  1  »  ) 

(C0L(1 ),A(4),X<6)*3  F<10.5)*X(6)»A)S 
GO  TO  A 1 S 
D 1 : END  fts 


LCD  (  I  )  * 


=  »tKSAT» 
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ATTACHMENT  C 
POTEET  Program 
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o  n  n  noo  n  r>  ooo  '■>"'>  r>  r>  o  n 


T?hl p  7. — POTEET  proavctr  listing 


C 


*»eeeaee4MM>a«oa»eoeefM>ee&aoae<M>oaaoo«0oaeeoaaaao4>a<MM>oa6<»aoaeaatt  POT 

POT 

POTEET  POT 

POT 

PUT 

THIS  PROGRAM  COMPUTES  DAILY  POTENTIAL  FVAPOTPANSPIRATIOM.  POT 

IN  INCHES  PER  DAY,  USING  A  METHOD  DEVELOPED  BY  C.  W.  ThORNTh m A  I TE . PO T 
PRIMARY  INPUT  IS  DAILY  MAXIMUM  AND  MINIMUM  AND  MONTHLY  AVERAGE  PQT 

TEMPERATURE  DATA  FROM  WEATHER  BUREAU  STATIONS.  POT 

POT 

POT 

COMMON/C?/ IRD. IPT . IPCH  POT 

REAL  MINT  (  1EIS0)  .MAXT  (  1BS0)  POT 

DIMENSION  AMT  (  1  2)  ,PE  (  1  85(1  )  POT 

I PD= 1  POT 

I  PC H  =  1  6  POT 

I P  T  =6  POT 

POT 

NSTAS=NUMBER  OF  WEATHER  BUREAU  STATIONS  FOR  WHICH  POT  FT  IS  COMPUTPOT 

POT 

READ  ( I  PD .11)  NSTAS  POT 

DO  S  I JKLMN=1 .NSTAS  POT 

POT 

AMT  =  AVFRAGF  MONTHLY  TFMPEP ATURF  «  IN  DEGREES  E  POT 

POT 

READ  (IR0.10)  < AMT < J> , J= 1 , 1?)  POT 

POT 

HTI=THORNTHWAITE  HEAT  INDFX  POT 

POT 

HTIrO.n  POT 

DO  3  J=l,12  POT 

IF  (AMT(J)-32.)  1.1.2  POT 

1  Hl=0.0  POT 

GO  TO  3  POT 

2  HI= (  ( AMT ( J) -32. ) /Q. ) «* 1 .SI  4  POT 

3  HT I =HT I wH I  POT 

POT 

A=THORNTHWA I TE • S  EXPONFNT  POT 

POT 

A= (6.75F-07* (HTI**3. )  > - (7.71F-0S* (HTI**2. >  >  ♦ ( 1 . 7RE-02*HT I )  »4 .BF-01 POT 
PI=3. 14160  POT 

POT 

ST AL AT=STAT ION  LATITUDE,  AS  A  DECIMAL  NUMBER.  POT 

POT 

READ  (TP0.9)  STALAT  POT 

POT 

AMP=AMPLITUOE  of  sine-wavf  variation  IN  DAYLIGHT  factor  POT 

POT 


AMP= ( 1 .86E-0S* (STAL  AT** 3. )  ) - (2.087E-03* (STAL  AT* *2. )  ) ♦ (8. SI 7E-0?*STPOT 


1 

2 

3 

4 

s 

6 

7 

8 
R 

1  0 
I  1 
1  2 
1  3 

1  4 

1  S 
16 
1  7 
18 
I  B 
20 
21 
22 

23 

24 

25 

26 

27 

28 
2B 

30 

31 

32 

33 

34 

35 

36 

37 

38 
3B 

40 

41 

42 

43 

44 

45 

46 

47 

48 
4  B 
50 


82 


1 


r%  o  o  •">  non 


Table  7. — PJTEET  program  listing — Continued 


1  ALAT) 

POT 

S] 

c 

POT 

5? 

c 

NYR=NUMREP  of  years  of  station  record  FOP 

THE 

GIVEN  STATION 

TO  HE  POT 

S3 

c 

ANAL V7F0 

POT 

S  4 

r 

POT 

S5 

RF  AO  < I RD  « 1 1 )  NYR 

POT 

SR 

r 

POT 

S  7 

r 

MOrTOTAL  NUMBER  OF  DAYS  IN  PERIOD 

TO  Hf  PROCESSED. 

POT 

SH 

c 

POT 

SR 

r 

OSE  =  OA  YS  SINCE  SPRING  EQUINOX  TO 

HFG INNING 

OF 

RECORD  TO  HE 

ANAL  YZEPOT 

GO 

c 

(-80  OR-81  FOR  JANUARY  1) 

POT 

G  1 

c 

POT 

G? 

READ  (IRD.R)  MO. OSE 

POT 

G  3 

r 

POT 

G4 

r 

I  YE AR  =  CALFNOAR  YEAR  FOR  WHICH  POTENTIAL  ft 

IS 

COMPUTE  1 

POT 

GS 

C 

POT 

G  G 

READ  (IRO.7)  IYFAR 

POT 

G7 

CALL  RFADSO (MAXT.MO. IRO) 

POT 

GR 

CALL  REAOSO(MINT.MO, IRO) 

POT 

GR 

pfsum=o. 

ROT 

70 

NO=wO- 00 

POT 

71 

00  4  1=1. NO 

POT 

7  ? 

K= I ♦ 3n 

POT 

73 

TSUM  =  MAXT <k ) .mint  <k I 

POT 

74 

TFMP= (TSUM/P.-3?. )  /  I  .8 

POT 

7  S 

IF  (TRMP.lT.O.)  TFmP=o. 

POT 

7  G 

c 

POT 

77 

r 

DLF  =  DA Y  LENGTH  FACTOR,  THF  RATIO 

OF  HOURS 

OF 

DAYLIGHT  TO  1? 

POT 

78 

r 

POT 

79 

0LF  =  l.n.((AMP-l.)*SlN(Pl»(I.f)SF)/l«3.)) 

POT 

80 

UPF=.0?1»< ( (10.*TFMP)/HTI)*»A) 

POT 

«1 

PF ( I ) =HPE  *0L  F 

POT 

8? 

& 

PFSUM  =  PFSl)M.PE  (  I  ) 

POT 

88 

s 

CONTINUE 

POT 

84 

WRITE  (IPCH.M  (PF ( I )  .  1  =  1 ,M0> 

POT 

85 

WRITE  (IPT.l?)  P  E  S  U  M 

POT 

8  G 

STOP 

POT 

87 

c 

POT 

88 

G 

FOPMAT  (  1  OF  7 . 4 , 3X . » PE •  ) 

POT 

89 

7 

FORMAT  (14) 

POT 

90 

8 

FORMAT  (I4.FS.0) 

POT 

91 

9 

FORMAT  (F8.7) 

POT 

9? 

1  o 

FORMAT  (1PFG.R) 

POT 

93 

1  1 

FORMAT  (I?) 

POT 

94 

1  ? 

FORMAT  (  1 X .  • Pfc SUM= •  ,F8  .  4 ) 

POT 

9  S 

END 

POT 

9G- 

8  i 


•ffiaiiriiiimtni-ifTr- 


r>  r>  n  ooooooo  ooor>  o  r>  n  o  o  o  o 


Table  7 .  — Par  EE?  vvogran  l  i  s  i  i  Kg  —  Continued 


SUBROUTINE  PEA0S0(S0. ICNT, IPO) 

REA 

1 

REA 

2 

REA 

3 

REAOSO  INPUTS  THF  NUMBER  OF  OATS  IN  EACH  MONTH  AND  MAX.  AND  MIN 

.REA 

4 

temperatures. 

REA 

5 

REA 

6 

REA 

7 

REA 

8 

DIMENSION  SO(l BSO) . MOA Y ( 1 ? ) .TEMP (I860) , IYR(IBSO) ,  MON (1850) , 

REA 

9 

1  I  PAY ( 1 850 ) 

REA 

10 

REA 

1  1 

MDAY(IO)  —  ARRAY  CONTAINING  THE  NUMRER  OF  DAYS  IN  EACH  CALENDAR 

PEA 

12 

MONTH. 

REA 

1  3 

REA 

14 

READ  ( T  RD  » 9 )  ( MDAY ( IQ) , I G= 1 , 1 2 > 

REA 

15 

1  =  0 

REA 

16 

REA 

17 

I YRD  —  CALENDAR  YEAR 

REA 

18 

REA 

19 

IMOND  —  CALENDAR  MONTH 

REA 

20 

REA 

21 

TEMP ( J )  --  TEMPERATURES  FOR  FIRST  10  DAYS  OF  MONTH. 

PEA 

22 

REA 

23 

1 

READ  ( I RD  «  7 )  I YRD . I  MONO  «  (T£MP{J)  . J  =  1 , J  0 ) 

REA 

24 

IF  (IMOND-2)  4,2.4 

REA 

25 

2 

I  X Y= I YRD/4 

REA 

26 

IF  (  ( I X  Y*4 ) - 1 YRD )  4,3.4 

REA 

27 

3 

MDAY (?) =29 

REA 

28 

4 

ISTOPM=MDAY (IMOND) 

REA 

29 

REA 

30 

TFMP(J)  —  TEMPERATURES  FOR  DAY  11  TO  END  OF  MONTH. 

REA 

31 

REA 

32 

READ  (IRD.B)  (TEMPI J) ,J=H .ISTOPM) 

REA 

33 

MDAY (?) =2B 

REA 

34 

DO  S  J=l, ISTOPM 

REA 

35 

1  =  1*1 

PEA 

36 

I YR ( I ) =1 YRD 

REA 

37 

MON ( I ) =IM0ND 

BFA 

38 

I 0  A  Y ( I ) =J 

REA 

39 

SO(I) =TEMP ( J) 

REA 

40 

s 

CONTINUE 

REA 

41 

IF  (I-ICNT)  1,6,6 

REA 

42 

* 

CONTINUE 

RFA 

43 

RETURN 

REA 

44 

REA 

45 

7 

FORMAT  (9X.2I2, IX, 10F6.2) 

REA 

46 

8 

FORMAT  (14X,10F6.2,/,14X,11F6.?) 

REA 

47 

Q 

FORMAT  (121?) 

REA 

48 

END 

REA 

49 

\ 

\ 


84 


'lumber 


Sequence  number  of  day  in  inpt 
master-data  set  where  interpo 


Table  9. 


*o  -» >i 


;■  .  ... 


R  I  --  kU' 

INTFppql  at  IONj  AND  avcijaaHn;  r  I  v 

program  Piv 

(FOR  MAIN ST FM|  R  I  V 

R  I  V 

p  I  V 

T  p  1  ^  PwOO^flM  IS  D F  S I GN F  o  To  P  R  0  V  I  O  F  T  mF  GRODNO- w  A  t  F  p  F  [  r  i  w  P  I  v 

SIMULATION  MOI1KL*  SURF  o^nf<  .  WITH  IG-DAY  AVFkATF  UIVF--FTAF.F  FATA  u  I  V 

FtfFFJy  10  r.AvS  FOP  P  P  F.  c  I  F  I  F  1  I  CORRF  G  POND  I  NT.  NOOF  LFvFIS  t*lO  Plvpc  PIV 
HU  FS.  P  I  v 

PRIMARY  INPUT  I  S  roPPFFPiiNr.^r,  -IvFR-STAGF  ANT  W[Kfw  mjlf  njTAPIv 
IN  f-0  a  y  I  AJC'.PF  MF  NTS  wh;th  tan  rF  ■jFj!)  Fi-Oy  F  I  T  H  F  P  maGuFTIC  ’apf  RIv 

op  niSK  FH.F’S  nw  FROM  CARDS.  NOF'F  LfVfLG  A  n  "i  TMF  I  w  A  PRROPR  I  ATF  P I  V  P I  V 

MILFS  ARF  PF  ATI  FROM  CAPPS.  Ptv 

THF  POOFFAM  I  N  T  F  P  p  0  L  A  T  F  S  Fpp  R  0  T  H  T  I  M  p  ANP  p  I  V  F  P  M  I  L  F  ^  AN1.1  F  I  V 

COMPUTFS  10-OAY  A  V F  P  AGF  S  F(,r  IhF  FvTJRF  PFRJOO  OF  RFCron.  PIv 

THF  F  I  PS  T  RFCflFJf)  IN  T«F  OUTPUT  OAT  A  Sr  T  TFpLS  YOU  hD^  many  NOORIV 
LFVFLS  YOU  HAVF  RJvFR  S  T  A  Or  FOP.  how  many  (,«0'JD?  OF  1>,-DAY  ~Iv 

AVfcPAGFS  YOU  hSVE.  ANP  ThF  T  I  «F  INCRFmFNT.  IN  ['AYS.  -IV 

THF  DATA  APF  WRITTFN  ONTO  A  maGNFTIC  STORAGF  I/FVICE  IN  Plv 

UNFORHATTFO.  V  AP  I  API.  F -IF  NO  I  h  pf'CORDS.  ulv 

p  I  V 

-  I  V 

•  [NPIJT  DATA  ooo  PIV 

-  T  V 

IMON  -  R  F  G I NN I  NO  month:  p  IV 

IDAY  -  1 0  A  Y  ♦  )0  :  day  OF  FIRST  1  ■-  IAY  AvFPAhF  I.IITRUT:  PIV 

IYFAR  -  BF ginning  yfar  r I v 

-I  v 

ICNT  -  NUMPFR  OF  DAYS  PFRIOTi  OF  R-  -O  0  O  v  F  -  S  .  P’v 


N  0  A  Y  S  -  I  IhF  InORFmFNT  r[V 

-  I  V 

NSTAGD  -  NIJMSt.o  OF  NU^f  Lt  Ul  F  RIv 

-  I  v 

u  -  ADRAY  HOLDING  NO  OF  LF'vFlS  R  I  v 

R  I  V 

OM  -  ARRAY  HOLDING  PIVFP  HUES  rf)-RF  Sr  I  NS  TO  NOOF  I.FVFI.s,  RJV 

R  I  V 

1 1 1  m  -  niimmfp  Of  corpf  spono  I  Nf-  rjvFr  mi;fn  and  -ivfr  stages  fop  fatriv 

h  TIME  STFP  in  THF  INPUT  "114  SFT.  SJV 


ISTART  -  SFOUFNCr  N'lMHFR  of  day  wf  L  ATIVF  f  O  INPUT  DA  T  «.  sf  T  k«i  RF  RJV 
COMPUTATION  to  RF-,IN.  P  I  A 

R  I  V 

IHFGN  -  SF  O'  IF  NT  F  NI'MRFR  OF  NAY  IN  INPUT  '  A  I  A  SfT  «RM'  I  N  T  F  P  -  o  1  A  T  I  R  I  , 


S  AD-A099  337 


UNCLASSIFIED 


GEOLOGICAL  SURVEY  BATON  ROUGE  LA  MATER  RESOURCES  DIV 
METHODS  ANO  APPLICATIONS  OF  DI6ITAL-M00EL  SIMULATION  OF 
MAY  80  A  H  LUDMI6»  J  E  TERRY 
USGS/MRD/MRI-81/037 


F/6  0/8 
THE  RED— ETC (U) 


d  d  d 


Table  9. — RIVCHANGE  program  listing— Conti nued 


r 

c 

c 

c 

c 

c 

c 

c 

r 

c 

r 

c 

c 

c 

c 

c 

r 


r 

c 

c 

c 

c 

r 

r 

r 

c 

c 

c 

c 

r 

r 

c 


TO  OEGIN. 


TO  FND. 


OUMMY  -  DATE  ON  INPUT  DATA  SET. 

DAY  -  SFOUFNCE  NUMflFR  ON  INPUT  nATA  SF  T  (MllLTIRLFS  OF  F I VF >  . 
6mm  -  ARP  A  Y  HOLDING  R  T  V  F  R  MILES  ON  INPUT  DATA  SF  T  * 


INPUT  DATA  SFT. 

***  OUTPUT  DATA  *** 

NSTA6F  -  Niimrfp  OF  NOOF  LFvFLS. 

NOSFT  -  N'JMPFR  OF  GROUPS  OF  10-OAY  AVERAGES  (RECORDS)  YOU 
HAVF  IN  OUTPUT  DATA  SFT. 

NDAYS  -  T I MF  INCREMENT 

IRON. IOAY. I YFAR  -  MONTH,  DAY,  AND  YEAR.  FACh  OUTPUT  RFfOPD  IS 
DATFO. 

J  -  SFOUENCF  NUMRFP  OF  FACh  RECORD  IN  OUTPUT  DATA  SET. 

Iu  -  ARRAY  HOLDING  NOOF  LFVFLS  IN  OUTPUT  DATA  SET. 


IN  OUTPUT  DATA  SET. 

DIMENSION  JJ{?0n),RM(pnn),H(?O0).GH(S0.Pno),GM(S(1),FI.  (Sn,?) 
DIMENSION  JO  ay  (]  ?)  ,GMM  (1  SO  )  ,FEL  (iso  .?)  .oiimmy  (  ?(i) 

DATA  JDAY/31  ,?P  ,  31  .  TO.  11  *10,  il  ,  31 , 30.D1  »3f,  ,  11  / 

DATA  GH/ 1 0000*9  ,/ 

I  RD  =  S 
I  TP=  1  0 
IDA  =  P 
IPT  =  S 

READ  t  T  RO  »  1  ft )  IRON,  IDAY  .  I  Ye:  AR 

IF  (MOD( IYFAR.4) .FQ.O)  jnAY(?)=?P 

RF  An  ( I RD ,1Q)  ICNT, NDAYS, NSTAGE 

RF  AD  (  T  BD  »?()  >  (  1  J  (  I)  ,RM  (  I  )  ,  1  =  1  ,NST  AGP  ) 

RF AD  (T»D*?4>  IUM 
READ  (IRD.P4)  ISTART 
RF  AD  ( I RD • ?4  >  IPEGn.IEND 


PI  V 

G) 

PI  V 

s? 

Civ 

S3 

RIV 

S4 

R  I  V 

SS 

p  I  V 

SF> 

PI  V 

S  7 

PI  V 

SS 

Piv 

SR 

PI  V 

so 

PIV 

si 

R  I  V 

S? 

Piv 

S3 

«lv 

S4 

Piv 

S  S 

PIV 

SS 

R  1  V 

b  7 

PIV 

Sft 

PIV 

SO 

R  1  V 

70 

PIV 

71 

P  I  V 

7? 

PIV 

73 

PIV 

74 

PIV 

7  S 

PIV 

7  S 

RIV 

77 

PIV 

7ft 

PIV 

70 

PIV 

SO 

PIV 

«1 

PIV 

ft? 

PIV 

ft  3 

PIV 

ft* 

PIV 

ftS 

Piv 

PS 

Piv 

ft  7 

PIV 

flft 

PIV 

ft  P 

PIV 

O0 

PIV 

PI 

P  I  V 

P? 

p  I  V 

P  3 

RIV 

P* 

RIV 

PS 

R  I  V 

PS 

PIV 

P  7 

P  I  V 

Qft 

R  I  V 

OP 

P  I  V 

1  00 

89 


Table  9. —  RIVCHANGE  program  listing — Continued 


1  READ  (ITP.21)  DUMMY 
READ  ( I  TP  »  25 )  DAY 
1 1 OAY  =  DAY 

READ  (IT P,?A)  (GMM  <I>  .EEL  ( I « 1 ) » 1=1 » IUM) 

IE  (I10AY.NE.I8EGN)  GO  TO  1 

JZ  =  0 

IZ  =  0 

DO  2  J=1.IUM 

IF  (GMM(J) .GE.RM(l) .AND.IZ.EG.O)  IZ=J-1 
IE  (GMM(j) .GE.RM(NSTAGE) . AND. JZ.EQ.O)  JZ=J 
?  CONTINUE 
NUM=0 

DO  3  I=IZ.J7 

NUM=NUM*1 

GM(NUMJ=GMM(I) 

3  El ( NUM , 1 ) =EEL ( I  »  1 ) 

CDAYS=ISTART-I 1DAY 
IXIN=ISTART 

IS  = ISTART/10 

4  READ  (ITP.21)  DUMMY 
WRITE(IPT,201)DUMMY 
READ  ( I  TP  »  25 )  DAY 
WRITE ( IPT  ,  30 ) DA Y 

RE  AO  (ITP.27)  (EEL ( I .2) » 1=1 » IUM) 

NZ=IZ 

DO  5  J=1,NUM 
EL(J.2)=EEL(NZ,2) 

5  NZ=NZ*1 

WRITE  ( IPT  »  30)  (EL ( I ,2) » 1=1 »NUM) 

I  PDA Y=DA Y 
IX=I2DAY-I1DAY 
I D 1 = I 1 0 A Y* 1 

IE  (IOl.LT.ISTART)  ID1=ISTART 
DO  6  J= ID  1 , I 20A Y 
JJ  =  J- 1 1  DAY 

JTH=(J-ISTART)/10*IS 
DO  6  1=1 .NUM 

DEL*(EL(I.2)-EL(I,1> )/IX 

6  GH  ( I  «  JTH)  =  EL  (1*1)  *DEL*JJ*GH( I , JTH) 

DO  7  1=1. NUM 

7  EL ( I . 1 ) =EL ( I *?) 

I10AY=I2DAY 

IE  (DAY-IEND)  4,8.8 

8  DO  9  1=1.50 
DO  9  J=1 ,200 

9  GH ( I , J ) =GH ( I , J ) / 1 0 • 

NQSET=  TCNT /NDAYS 

WRITE  (IDA)  NSTAGE, NOSET, NDAYS 
WRITE  (IPT, 28)  NSTAGE. NOSET, NOAYS 
ISTATH=ISTAPT/10 


RIV  101 
RIV  102 
RIV  103 
RIV  104 
RIV  105 
RIV  106 
RIV  107 
RIV  108 
RIV  109 
RIV  110 
RIV  111 
RIV  11? 
RIV  113 
RIV  114 
RIV  115 
RIV  116 
RIV  117 
RIV  118 
RIV  119 
RIV  120 
RIV  121 
RIV  122 
RIV  123 
RIV  124 
RIV  125 
RIV  126 
RIV  127 
RIV  128 
RIV  129 
RIV  130 
RIV  131 
RIV  132 
RIV  133 
RIV  134 
RIV  135 
RIV  136 
PIV  137 
RIV  138 
RIV  139 
PIV  140 
RIV  141 
RIV  142 
RIV  143 
RIV  144 
PIV  145 
RIV  146 
RIV  147 
RIV  148 
RIV  149 
RIV  150 


Table  9 .—RIVCHANGE  oroavan  listing  —  Continued 


r 


00  17  J= 1ST ATH # JTH 
00  13  1=1 .NSTAGE 
00  10  K=1 »NUM 

IF  <RM(I> .GE.GM(K) )  GO  TO  10 
GO  TO  1 1 

10  CONTINUE 
K  =  NUM 

GO  TO  12 

11  IF  (K.FQ.l)  GO  TO  12 
KS=K-1 

H ( I ) = ( RM ( I ) -GM ( KS  > )/(GH(KS*l)-GM(KS) ) * ( GH < KS* 1 . J> 
1  J) 

GO  TO  13 
1?  H ( I ) =GH (K  ,  J ) 

13  CONTINUE 

IOA Vs IOAY  *10 

IF  (IOAY.LE.JOAY(IMON)  )  GO  TO  16 
IDAY*I0AY-JDAY ( IMONI 
I MON= I MON* 1 

IF  (IM0N.LE.12)  GO  TO  16 
I MON= 1 

iyear=iyearu 

IF  (M00< IYEAR.A) )  15.1*. 15 
1*  JOAY (2) =29 
GO  TO  16 

15  JOAY (2) =20 

16  WRITE  ( IOA )  I MON .IOAY.IYEAR 
WRITE  ( IPT.22)  IMON, IOAY.IYEAR 
WPITE  (IDA)  J 

WRTTE  (IPT.29I  J 

WRITE  (IOA)  (IJ(I) .M(I) ,I=1.NSTAGE) 

WRITE  (IPT.23)  (IJ(I) .MCI) ,1*1 .NSTAGE ) 

17  CONTINUE 


STOP 

18 

FORMAT 

(312) 

19 

FORMAT 

(1615) 

20 

format 

<8(I*»F6.1) ) 

21 

format 

( 20A* ) 

2? 

FORMAT 

(1X,I2,»/».I2.‘/».I2I 

23 

FORMAT 

(TUX. I*. IX. F6. 2)  ) 

2* 

format 

(215) 

25 

FORMAT 

(F10.3) 

26 

FORMAT 

(8F10.3) 

27 

format 

(♦(10X.F10.3) ) 

28 

FORMAT 

(IX. 515) 

29 

format 

END 

(IX. 15) 

RIV  151 
RIV  152 
RIV  153 
RIV  15* 
RIV  155 
RIV  156 
RIV  157 
RIV  158 
RIV  159 
RIV  160 

-GH(KS.J) ) ♦GH(KS,RIV  161 
RIV  162 
RIV  163 
RIV  16* 
RIV  165 
RIV  166 
RIV  167 
RIV  168 
RIV  169 
RIV  170 
RIV  171 
RIV  172 
RIV  173 
RIV  17* 
RIV  175 
RIV  176 
RIV  177 
RIV  178 
RIV  179 
RIV  180 
RIV  181 
RIV  18? 
RIV  183 
RIV  18* 
RIV  185 
RIV  186 
RIV  187 
RIV  188 
PIV  189 
RIV  190 
RIV  191 
RIV  192 
RIV  193 
RIV  19* 
RIV  195 
RIV  196 
RIV  197 
RIV  198 


91 


Depends  1-75  15F5.2  GH(I,J),  Array  holding  input  10-day-average 
on  J=1,NDATA,  stream  stages  for  each  gage  on  the 
NDATA  1=1 ,NGAGES  stream. 


Table  11 s,  'HAX  IF  nr,  ’;)r,rl 


1M 

T  R  1  e>ChangE  --  TPI 

IMF  PPOL  AT  I  ON  PROGRAM  TUI 

(FOP  TPlHl'TAPY  ST»H«S1  TRI 

TUI 


TrI 

THIS  PROGRAM  OUOVIDFS  1 P-DA Y- AVE PAGE  STOE  am-STAOE  •'lATA  TP  I 

OOPPE  SPONO INO  TO  SPFOIEIC  NODE  LEVELS  For  T.-IhuIary  SIrEA-'S  T  >-  I 

IN  The  awe  A.  T  R I 

1 O-DAV-AVtPAGE  OATrt  EOS  OlEEfRrNT  Eaiqw,  SITES  n.  MCH  T“I 

STREAM  ARE  RE  All  .  I  NT  f  DPOI  A  T  I  ON  f  OP  OISTANOE  IS  RE  RE  U  m  f  f  •  T  W  1 

IN  OROFW  TO  DETERMINE  STAGES  EOP  PAuT|CU|AR  NO^ES.  TRI 

OUTPUT  IS  WRITTEN  ON  MARNE  TIT  <'[s«  EUe  TO  me  REEERENOEU  RY  TRI 
MOREL  I  NO  PROGRAM.  TRI 

TRI 


iaeoeo44aQani}»u«oeDOi.'i]oe>>^ 


4  44  4D4444  44444444444  >4  4  4  44444  J  4  444  J4  44  44444444  >44444  ThI 


TRI 

«»»  input  DATA  *«<*  T  R 1 

TRI 

DIIM  -  TITLE  he  a  o  i  Nfi  E  or  prime*'  output.  TRI 

TR[ 

NS T RMS  -  numrer  OP  TRT«UTary  STREAMS.  TPJ 

TRI 

NOATA  -  NUMRfw  OF  lrt-»1AY  AVERAGES  RE  I  JO  ENTERED  f  OR  E  1  f  H  STREAM.  T  R  1 

TRI 

NAVE  -  NUMBER  OF  m-DAY  AVERAGE  RECORDS  To  RE  ]N  OUTPUT.  TrJ 

T«I 

NS  T  0  T  -  TOTAL  NUMBER  OF  NOOES  TO  me  ASSTGnE”  A  STREAM  STAOE.  TRI 

TR  1 

—  REPETITIVE  FOR  EACH  STREAM  TR I 

TRI 

NOAO,F  -  number  OF  GAGES  ON  STREAM.  TRI 

TRI 

NSTAGE  -  N'lMREP  OF  NOOES  APPL  I C  AHI  E  To  T«;s  STREAM.  T»I 

TRI 

GEl  -  ELEVATION  OF  GAGF.  TRI 

TRI 

r,M  _  STPF  AM  MU  F  OF  GAGE  .  TRI 

TRI 

GH  -  ARRAY  HOLOINO  1  0 -(1 A  Y  -  A  Vf  R  A  OE  STR1am  STAGES.  T  ‘M 

TrI 

IJ  -  ARRAY  HOLDING  NODE  LEVELS.  T«I 

TRI 

RM  -  ARRAY  HOLDING  STREAM  Mil  f  S  COPRE  SE'omo  I  NO-  To  note  LEVEIS.  TR] 

TRI 

*«*  OUTPUT  DATA  ***  TRI 


Table  11.  —  TRIBCHANut'  vroaran  listing — Continued 


r 

OHM  -  TITLE  HEADING  (IN  PRINT  ON|_  Y  )  . 

T  P  I 

SI 

r 

T«I 

S? 

c 

NS  TOT  -  TOTAL  NUMHfP  OF  NODES. 

T  P  I 

S3 

c 

T  P  I 

S4 

c 

NAVE  -  NUMRfP  OF  1  0 -0  A  Y  -  A  VE  P  A  GF  RECORDS  IN  OUTPUT  (PRINT  ONIY). 

T  P  I 

GS 

r 

T  P  1 

SG 

c 

IJ  -  APPAY  HOLDING  ALL  NOPE  S . 

T  P  I 

S  7 

r 

7  P  I 

SR 

c 

k  -  SEQUENCE  NUMHFP  FOR  Each  RECORD  wpITTF1". 

T  P  1 

SR 

c 

I  P  I 

GO 

c 

H  -  APPAY  HOLDING  STREAM-STAGE  VALUES  CORRESPONDING  Ti  SPECIFIFD 

T  p  I 

G  1 

c 

NODE  LFVELS  AND  STPFAM  milFS. 

T  P  I 

G  A 

r 

TP  1 

G3 

r 

°  T  P  I 

G4 

r 

T  P  I 

GS 

DIMENSION  GH  <  7. 300) . I J  <  300 ) .R«  <  300 ) ,n  <  ?o 1 . 1 R?) .GM ( 7) .GE  l  <  7  > 

T  P  I 

GG 

DIMENSION  DUM ( ?  0 ) 

T  R  1 

R7 

OATA  h/3^s«P*I)./ 

T  P  I 

GR 

OATA  JPO/S/.IDA/i?/.  IPT/*./ 

T  P  1 

G  9 

READ  (IPO.O)  Ohm 

T  P  I 

70 

READ  (TRD.10)  MSTRMS.NOATA 

T  P  I 

7  1 

RFAD  < I RO  « 1 0  >  NAVF.NSTOT 

T  P  I 

7? 

I  T  T  =  0 

T  P  1 

7  I 

DO  7  I S= 1 , mSTpmS 

T  P  I 

74 

RFAD  ( I R  0 . 1 0 )  ngagfs.nstage 

T  p  I 

7  S 

DO  1  I  =  1 <  USAGES 

T  P  I 

7R 

READ  (IRO.ll)  GFL(I).GM(I) 

Txl 

77 

RfAD  (TRD.ll)  (GH (  I  , J)  , J=1  .MDATA ) 

T  p  1 

7  R 

1 

CONTINUE 

7  P  I 

7  9 

I  T= I TT  »  1 

T  P  I 

RO 

ITT=ITT*NSTAGE 

T  P  1 

R1 

READ  ( I RD . 1 ? )  (  I J ( I )  ,om(I) ,I  =  IT.ITT) 

T  p  I 

R  ? 

DO  R  j= 1 . NA VE 

TP  I 

a  3 

DO  S  I=IT,ITT 

T  p  1 

R4 

DO  ?  k= 1 . MG AGE S 

T  P  I 

OS 

IF  (PM(T ) .GF.GM(K) )  GO  TO  ? 

T  P  1 

RG 

GO  TO  3 

TP  I 

R  7 

? 

CONT  I  NIIF 

T  P  I 

MR 

KrNGAGFS 

T  P  I 

RR 

GO  TO  A 

T  P  I 

Q0 

IE  (K .FO. 1 )  GO  TO  S 

T  P  I 

R  1 

K  S  =  K  —  1 

T  P  I 

R? 

H(I.J)  =  (WM(I)-GM(K’S))/(GM(KS*l)-(iM(KS))®(GH(KS.l.J)-Gwi(KS,,l).GEL(('TPl 

Q3 

1S*1  l-GFL  (KS)  )  *GF  L  (KS)  *GH(KS»J) 

TP  I 

P  4 

GO  TO  S 

T  P  1 

PS 

4 

IE  (RM( J ) ,GT ,GM(K ) )  GO  TO  S 

T  P  I 

RG 

H (  I  . J) =GH (K. J> ♦ G  E  L ( K ) 

T  P  I 

R  7 

c, 

CONTINUE 

T  P  I 

OP 

G 

CONT INHF 

T  P  I 

RR 

7 

CONT  INIlF 

T  P  I 

1  no 
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Table  11.—  TRIBCHANGE  rroarrjm  1  t in.; — Continued 


WRITE 

1IPT. 13)  DUM 

TRI 

101 

WRITE 

(IPT.14)  NSTOT*NAVE 

TRI 

102 

WRITE 

(IDA)  NSTOT 

TRI 

103 

WRITE 

(IDA)  (IJ(I) .1=1, NSTOT) 

TRI 

104 

WRITE 

(IPT.U)  NSTOT 

TRI 

105 

WRITE 

(IPT.1S)  (Util  .1  =  1. NSTOT) 

TRI 

106 

K  =  9 

TRI 

107 

DO  9  J 

=  1  .NAVE 

TRI 

108 

WRITE 

(IPT.U)  K 

TRI 

109 

WRITE 

(IDA)  K 

TRI 

110 

WRITE 

(IPT.U)  (IJII)  ,H(I,J)  ,1  =  1, NSTOT) 

TRI 

111 

WRITE 

(IDA)  (IJ(I) ,H(I,J) ,1=1, NSTOT) 

TRI 

112 

K=K  ♦  1 

TRI 

113 

8 

CONTINUE 

TRI 

114 

STOP 

TRI 

1  IS 

TRI 

116 

9 

FORMAT 

(20A4) 

TRI 

117 

10 

FORMAT 

(215) 

TRI 

118 

1 1 

FORMAT 

(15FS.2) 

TRI 

119 

1? 

format 

(8(U«F6.1)  ) 

TRI 

120 

13 

FORMAT 

( 1X.20AA) 

TRI 

121 

1* 

FORMAT 

(IX, 215) 

TRI 

122 

IS 

FORMAT 

(20 ( 1 X . U  )  ) 

TRI 

123 

If 

FORMAT 

(10(1X,U,1X,F6.1)  ) 

TRI 

124 

END 

TRI 

125 
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ATTACHMENT  F 
DELETDELH  Program 


99 


Program 


Table  13. —  DELETDELH  prog  ran  listing 


OFLET  :  /*  «  «  s  »  s  «  «  «  »  «  »«i!  e  #  » .  k  a  9  c-  iMf  » if  » i)  ito  ■>  iKiv  »  (to  «■  *  it  c  *  s  it  it 

DFLTAET  /  I'lELTAH 

THIS  PROGRAM  COMPUTES  nt'LTA  FT  /  ■'if- 1  T  A  h  USING 
THE  R  I  DPL F  F 1 1 N  C  T  I  0 N 6 1  . 


«»(HHf»«HHt»iHHHt<itti>ittHHtHHiit-Jit»v«iHHHf«eit«i;»  if  if  it  it  if  if  if  if  if 

■>/ 

PROCEDURE  OPT  TONS  MAIN)  : 

OFCLAPE  TT|?:S,30).FwHO(Tn)  Of  rno)  ,»'FLl.'iO  ; 

ON  FNOFTl  F  (SYSIN)  00  TO  F N 0 i ; 

> « 

RFAO  VALUES  of  ET/SAT.HYO.CONO.  fop  DFPTmc;  OF  ]  TO  7  a 
FEET  ABOVE  The  wATER  TARI.F  F  0°  E  OOP  RANGED  IN  VFPTICAL 
HYOF.AULIC  CONDUCTIVITY. 

'  / 

GFT  FILE  (SYSIN)  FO  I  7  (  ET  )  (  COI  (  15.10  F  (  7  .  h  )  ,  x  (  1  (, )  )  : 

/* 

READ  DATA  FOP  INDIVIDUAL  OBSERVATION  wf  LL  S  -  ID.  NUMRpP. 

VERTICAL  HYDRAULIC  CONDUCTIVITY  FOR  maTFR]AL  from  (.AND  SURFACE 
TO  THE  WATER  TABLE  AND  FROM  THE  wATFR  7  A  AL  F  70  TOR  OF  The 
AOUIFFR,  THICKNESS  FROM  l  AND  SURFACE  TO  TOP  of  AOuIffP. 
and  AVERAGE  OFPTH  TO  water. 

--■/ 

INI:  GET  FILF(SYSIN)  EDIT  (  WFl  L  NO  ,  hCU  .  HCL  .  T  h  I  r*  ,  ot  *  )  (  col  <  ]  )  ,  a  <  a  )  .  y  (  <■ )  . 
A  f  ( 1  n  )  )  : 

IF  H  C  U  < • 0  A  THFN  005 

IF  HCU<.nOA  THEN  IF  XP=  ?  5 
FLSF  IFXPS35 
ENn : 

flse  no: 

IF  H  C  U  < . A  THFN  IEXP=A5 
FLSF  IEXPrSJ 
FND  5 

x .  F .  f  =  n .  : 

no  i  =  i  to  3n : 

E,F  =  o.  ; 

x  =  X  ♦  1  .  ; 

Y  =  X  5 
j=T  : 

A  1  :  IF  J>10  THEN  ETO=0.  : 

FLSF  ET0  =  HCU*ET ( IF xp, j)  : 

IF  THICK>Y  then  no; 

FLOW  =  HCL-»  (Y-X  )  /  (  Thick -Y)  5 
IF  FTO>FLOW  THFN  DO! 

y  =  y  ♦  l  . ; 

J  =  J* 1  : 

F  =  F TO S 
F  =F  L  Ow : 

oo  to  ai: 
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Table  13.  —  PFLETDi'LU  frog ram  listing — Continued 


ENOS 

ELSE  OOt 

G=E-F-ETO*ELOWS 

IF  G> 0 •  THFN  GWETO ( I > =F- (E-ETO) * (E-E ) /G8 
ELSE  GWETO ( I ) =0 • 5 
GO  TO  APS 
ENOS 

ENOS 

ELSE  DOS 

IF  E>. 0000005  THEN  DOS 
EETO=E-ETOS 

EHYT=ETO>HCL*EETO* (Y-THICK) S 

OY= ( SORT ( EHYT**?>4 . *HCL* (THICK-X)*EETO)-FHYT) 
/ <?.»EET0)  $ 

GWETO ( I ) =ETO*EFTO* ( Y-THIC* ♦DY) S 
ENOS 

ELSE  GWETO ( I ) =ETOS 
ENOS 

A?:  IF  GWETO (  I  )  > . 00P?2  THEN  GWETO < I) = . 00R22 5 

ENOS 
J  =  DTW*  «5  S 
00  1  =  1  TO  30  S 

IF  J<=1  THEN  GWETOJ=GWETO(l ) 5 

IE  J>30  THEN  GWETOJ=0. S 

IF  J>1«.J<  =  30  THEN  GWETOJ=GWETO < J ) S 

JF  I-=J  THEN  DET < I ) = (GWETO ( 1 ) -GWETOJJ / ( J-r ) S 

FLSF  DOS 

IF  j=l  THEN  OET(I)=GWETO(1)-GWETO(?) 5 

IF  J> 1 F J<30  THEN  OET(I)=(GWETO(J-l)-GWETO(J*l))/?.S 
IF  J  =  30  THEN  OFT  ( I )  =GWETO  (  ?9  )  -GWETO  (  30  )  5 
ENOS 

ENOS 

PUT  FILF (SYSPPINT)  ED  I T ( WELLNO * » HCU=  **HCU*'HCL=  »»HCL* 

•  AVE.  DTW  =  •  ,DTW* 'THICKNESS  =  ' *  THICK  * 

•  OTW(FT)  ET(ET/OAY)  DET/DHU/DAY)  *, 

<!,GWFTO(I) .OET(I)  DO  1=1  TO  30)) 

(PAGE • A (4)  *4  (  S*<IP  U  )  *  A  *  E  (10*5)  )  *  SKIP ( ?)  *A, 

30  ( COL ( 5 )  *F (?) *X  (3)  *F (10*7)  »  X (5) »F (10*7) ) )  S 
GO  TO  INIS 
FNOl (END  0ELET8 
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HCU=  0 

.02000 

HCL=  0 

.00500 

AVE.  DTW 

=  3.00000 

THICKNESS 

=  14.00000 

OTW (FT) 

ET (FT/OAY) 

OET/DH ( 1/DAY) 

1 

0.0021079 

0.0002980 

2 

0.0017863 

0.0002745 

3 

0.0015119 

0.0002880 

4 

0.0012104 

0.0003014 

5 

0.0009873 

0.0002623 

6 

0.0007551 

0.0002523 

7 

0.0005939 

0.0002295 

8 

0.0004560 

0.0002112 

9 

0.0003467 

0.0001942 

10 

0.0002647 

0.0001782 

11 

0.0002044 

0.0001634 

12 

0.0001601 

0.0001502 

13 

0.0001266 

0.0001385 

14 

0.0001020 

0.0001282 

15 

0.0000832 

0.0001191 

16 

0.0000686 

0.0001110 

17 

0.0000572 

0.0001039 

18 

0.0000482 

0.0000976 

19 

0.0000410 

0.0000919 

20 

0.0000352 

0.0000869 

21 

0.0000304 

0.0000823 

22 

0.0000264 

0.0000782 

23 

0.0000232 

0.0000744 

24 

0.0000204 

0.0000710 

25 

0.0000180 

0.0000679 

26 

0.0000160 

0.0000650 

27 

0.0000144 

0.0000624 

28 

0.0000128 

0.0000600 

29 

0.0000116 

0.0000577 

30 

0.0000104 

0.0000556 

Figure  29.— Example  of  output  from  DELETDELH  program. 
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